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Transonic Wind Tunnel Test Sections’) 
By PAUL F, MAEDER and Sines D. Woop, Providence, R. I., USA?) 


1. Introduction 


Since it is generally impossible to predict the performance of an aircraft 
by purely theoretical means, the Aerodynamics Engineer is accustomed to 
relying on experimental data obtained in wind tunnels. Theoretical reasoning 

_ is then used to establish the laws of similarity between the flow about the mudel 
in the wind tunnel and that about the actual aircraft. : 

As long as the desired flight speeds are well below sonic speed, no particular 

_ difficulties are encountered by testing the models either in a closed channel or 
in a free jet of air. The significant similarity parameter then is the Reynolds py 
number, which can be adjusted by varying either the speed or the dont of 

the airstream. 

However, at speeds close to or exceeding the speed of sound, a second 
similarity parameter, the Mach number, must be considered. Tests now must . 
be performed at particular Mach and Reynolds numbers in order to obtain 
complete similarity. : 

With the present day emphasis on such high speed flight the wind tunnel, “ 

_if it is to be of practical use, certainly must furnish correct aerodynamic data ie 
_in the transonic speed range. It is particularly desirable to have such reliable 
experimental results in this range since, due to the nonlinearity of the theore- 
tical problem, the design of transonic aircraft based on theoretical work alone 
is at this time practically impossible and, at best, will be extremely difficult | 
in the future. 

The principal difference between the flow about an aircraft flying in the | 
atmosphere and the flow about its model in the wind tunnel, if the aerodynamic ~ 
similarity parameters are the same, is caused by the finite extension of the 
tunnel airstream. It can be shown that in a conventional wind tunnel this dif- 
ference increases as the speed of sound is approached. Moreover, if a test 

- section with solid walls is used, at a certain subsonic Mach number the model. ~ 
will cause the same effect as the throat of a Laval nozzle. The speed of the i, 
‘upstream flow may be increased no further, and the tunnel is then said to be 
choked. Supersonic speeds in such a test section can only be attained by inser- 
ting a Laval nozzle upstream of the model which has a throat considerably 
1) This work was sponsored by the Office of Scientific Research of the Air Research and Re 
’ Development Command of the United States Air Force. iy 


_?) Brown University. 
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smaller than the smallest cross-section at the location of the model. Thus for 
an appreciable section of the transonic speed range no testing is possible in a 
closed tunnel. Since the range of Mach numbers in which tests are impossible 
decreases with increasing tunnel size as 1/(1 — M3.) for a given model, where 
M,, is the undisturbed Mach number, huge tunnels with tremendous power 
requirements would be necessary to obtain useful transonic data, if there were 
no special transonic test sections available to reduce the wall interference and 
eliminate choking. 

In the following paragraphs the principles of various transonic test sections 
will be discussed. 


2. The Influence of the Walls on Data Obtained in a Wind Tunnel 


Consider the streamline pattern of the flow about some body in free flight. 
It consists of a parallel flow of infinite extension into which the body is placed. ~ 
The body causes the streamlines of the parallel flow to be deflected in such a 
way that the mass flow inside a stream tube remains constant, and the centri- — 
fugal forces caused by the streamline curvature are in equilibrium with the 
pressure forces. 

Assume now that the infinite parallel flow is replaced by a finite stream 
surrounded by some boundary, as will be the case in a wind tunnel. Then, in 
general, the streamline pattern about the body will not stay the same. Ifasolid — 
straight wall is used as the wind tunnel boundary, the streamlines forming the 
flow about the body are squeezed together more than they would be in free 
flight. This wall interference introduces an error and leads in some instances 
to the phenomenon of choking as already mentioned. The calculation of such 
interference effects at high subsonic speeds is discussed in [1], [2], [3]). 

If the tunnel is of the open jet type, i. e. the boundary is essentially one of 
air at rest, which is unable to support pressure differences, then the curvature 
of the outside streamlines becomes greater than for an infinite free stream in 
order to balance the forces caused by the body, since there is no longer any 
outside flow to resist deformation. A flow pattern will be obtained in which 
the streamlines are further apart than for free flight. The interference for such 
a high-speed open jet has also been determined (see, for example [14]). 

Thus, if interference is to be reduced or eliminated, some more complicated 
boundary which simulates the behavior of the free flight streamlines point for | 
point must be found. However, a perusal of the situation shows that this ideal 
boundary may not always be required in order to obtain essentially inter- 
ference-free data. 

Consider the disturbance induced by the boundary at a point on the body 
to be tested. If the speed at this point is subsonic, then it will receive inter- 


3) Numbers in brackets refer to References, page 209. 
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ference signals from a considerable portion, if not ail, of the boundary. Thus 
the total interference effect for this point is an average of the signals produced 
by the various portions of the boundary, and therefore the boundary has only 
to simulate the free flight streamline in the mean. 

If the speed at a point on the body is supersonic and the point cannot be 
reached from the wall along a characteristic line (i. e., the characteristic line 
extends only to the sonic line) the same averaging process takes place along 
the sonic line. The only time such averaging does not take place is when the 
speed at a point on the body is supersonic and it can be reached from the wall 
along a characteristic line. The wall then must simulate the free flight stream- 
line exactly at the origin of the characteristic line in order to yield zero inter- 
ference on the body. Even in this latter case it is conceivable that the point 
on the body is influenced by a great many characteristics originating at different 
parts of the boundary in such a way that no noticeable interference is produced. 
Whether or not this is true for a specific boundary will be discussed later in this 
paper. For the time being this case shall be excluded from the discussion. 

If the body is not too large, it is usually sufficient to assume that the inter- 
ference effects at any point on its surface are equal to the sum of the effects 
at some other point plus their linear variations in the x-direction between the 
two points. That is, the interference is approximated by the first two terms of its 
Taylor Series expansion. Letting « and v represent the velocity components 
parallel and perpendicular to the wall, respectively, this result may be expres- 
sed as: 


OUin 


(Uin)ent = (Uin)e-o + / ( ae Sei te parte wd Peers ( sit) 


This is the classical method of investigation and is based on determining, 
at only one point on the body, four types of interference; blockage interference, 
U;,; angle of incidence interference, v;,; buoyancy interference, 0u,,/0x; and 
streamline curvature interference, 0v;,,/0x. Experience has shown that if one 
accounts properly for these four types of interference, sufficient agreement 
between wind tunnel data and free flight data is obtained. 

Blockage interference is caused by the effect of the boundary on the 
distance between streamlines. In subsonic flow one can usually’compensate for 
this effect by applying a correction factor. The body behaves as if it were 
tested at a higher or lower speed than that measured in the tunnel, depending 
on the type of boundary used. In transonic flow, however, the tunnel may be 
either choked (solid wall), or a section of the desired speed range may be lost 
because of the change in character of the governing equations as the speed of 
sound is passed. This is, therefore, the type of interference which must receive 
the greatest amount of attention in the development of a boundary suitable for 
a transonic wind tunnel. 
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Angle of incidence interference is present if the slope of the streamlines is 
changed by the boundary due to its inability to counteract the lifting force 
of the body in the same manner as would the infinite airstream in free flight. 
For subsonic flow the well-known angle of attack correction can be applied to 
the test results to compensate for this type of interference. It is probably 
possible to apply a similar correction at transonic speeds. 

Buoyancy interference appears due to the fact that the boundary may 
cause a pressure gradient in the flow direction. This pressure gradient causes 
an additional force in the axial direction which is not present in the free flight 
condition. 

Streamline curvature interference will cause a symmetrical airfoil to be- 
have as if it were cambered, if tested at an angle of attack. That is, the effect 
of the walls is to change the streamline curvature to that which would be 
found for a cambered airfoil in free flight. 

Thus, except for points on the body which can be reached from the boun- 
dary directly along characteristic lines, it is sufficient to find a boundary which 
either eliminates or allows correction for these four types of interference at one 
point on the body throughout the transonic speed range. 

Since solid wall and open jet boundaries influence the streamline pattern 
in opposite manners, the solution obviously lies in a combination of the two. 
Early in the history of wind tunnels it was suggested [4], [5] that a tunnel with 
two solid and two open walls could be used to eliminate certain types of inter- 
ference. The next possibility to be investigated was the test section with one 
or two longitudinal slots [6], [7], [8]. PrstoLEsi [9], [10] and the NACA [11] 
went a step further by using walls with several longitudinal slots. Subsequent 
investigations with slotted test sections of various shapes [12]-(23] have shown 
that acceptable results can be obtained at transonic speeds. The idea of using 
perforations to produce the partly open effect [23], [31]-[37] evolved gradually 
from the use of porous walls [24]-[30] which, because of viscous effects, yield a 
homogeneous boundary condition for which the pressure drop across the wall 
is proportional to the mass passing through it. Although each of these transonic 
test sections utilizes the idea of a partially open boundary, a clear understand- 
ing of how the various boundaries are connected in their physical and mathe- 
matical aspects has been attained only gradually. 


3. The Boundary Condition Along Partly Open Walls 


A. General Remarks 


It will be helpful, before carrying out the mathematical analysis, to make 
a few general statements concerning partly open walls. As has already been 
mentioned, a satisfactory transonic test section must have some partly open 
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wall surface. The flow through the wall is controiled by varying the difference 
between the pressure P, of the plenum chamber on one side of the wall and 
the pressure P at a point in the test section which is far from the wall relative 
to the spacing of the wall openings. If AP = P — P, = 0 there is no flow 
through the wall, and the flow through the test section will have the undisturb- 
ed velocity U; if AP + 0, however, there will be a flow through the wall. 
The flow in the test section may then be found by superposition of a pertur- 
_bation flow onto the undisturbed flow. 

Without loss of generality this perturbation flow may be considered to be 
normal to the undisturbed flow at the openings. Thus the exact inhomogeneous 
boundary conditions become u’ = 0 at the openings and, since there is no flow 
normal to the solid portions of the walls, v’ = 0 elsewhere. These conditions, 
in turn, imply that the perturbation potential g’ is constant across the openings 
and the perturbation stream function yw’ is constant elsewhere. Since for 
longitudinal and swept slots these perturbation quantities vanish as x ap- 
proaches infinity, the conclusion is that g’ = 0 at the slots and y’ = 0 elsewhere. 
Now transverse slots are merely a special case of swept slots, and, as will be 
shown, the behavior of perforated walls is also similar to that of non-longitu- 
dinal slots. Thus the previous boundary conditions will be applicable to these 
cases also. 

At some distance far from the wall (relative to the spacing of the openings) 
the perturbation effects are averaged, yielding the velocities w’ and v’. Assuming 
an incompressible flow of density 9, Bernoulli’s equation between such a point 
and an opening (over which the pressure is P,) becomes 


P+ = [(U+w)?+ 0%) = P,+ 5 [U2 +04). (1) 


Neglecting the squares of the perturbation velocities, this equation simpli- 


fies to = 
AP=P—P,—=-e UU. (2) 


Thus the assumption of zero perturbation velocity in the axial direction at 
the openings requires that an axial perturbation velocity be introduced at 
some distance from the wall to account for the pressure difference across the 
wall. In actual practice the pressure difference across the wall for an empty 
test section will be zero, and the axial perturbation velocity will arise only 
when a body is placed in the test section. 

The original theoretical investigations to determine the interference in test 

sections with longitudinal slots (the so-called slotted wall test sections) were 
carried out by using at the walls the exact inhomogeneous boundary conditions 
just described [9]-{17]. These inhomogeneous boundary conditions led to a 
rather involved problem, since they caused basically simple flows, such as 
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two-dimensional or axially symmetrical flows, to become generally three- 
dimensional. In order to determine how slotted walls could be expected to reduce 
subsonic interference, either rash approximations had to be made, or extensive 
numerical work, usually requiring an electronic computer, had to be carried 
out to obtain the flow for even a few special cases. 

More recent investigations have shown, however, that it is sufficient in 
virtually all practical cases to replace the inhomogeneous boundary by a corres- 
ponding homogeneous mean boundary [18]-[21], [23] which, at some distance 
from the wall which is large compared to the distance between individual slots, 
will yield the same flow conditions as the inhomogeneous boundary (i. e., the 
uw’ and v’ already mentioned). If a great number of slots is employed, the 
introduction of such a mean condition is certainly permissible. 

The use of perforated walls has developed gradually from the concept of a 


porous wall. Since the flow through a porous wall is that of slow viscous motion, ~ 


a linear relationship between the pressure difference across the wall and the 


mean velocity component normal to the wall is obtained. Adopting the con- — 


vention that positive values of v’ represent an outflow through the wall, 
the approximation of equation (2) then leads to the mathematical condition 


Tae (3) 


where & is a positive constant of proportionality. 


If the diameter of the openings in the porous material is enlarged so that 


turbulent flow creates the pressure difference, the boundary condition assumes 
the mathematically inconvenient form 


w= +kd”? (4) 


the plus or minus sign depending on whether the flow through the wall is 


inward or outward respectively. 
In actual applications. the walls employed are relatively thin. The hole 


diameter is at least of the order of, if not considerably larger than, the wall 


thickness. The flow then resembles the orifice problem, which as is well known, 
can be explained by potential flow, neglecting viscous effects. Thus it was first 
assumed that the pressure in the hole cross-section would be proportional to 
the square of the normal velocity through it, leading to the same unfortunate 
nonlinear boundary condition as given by equation (4). It soon became obvious, 
however, that this is not the case, since the air will retain a considerable part, 
if not all, of its undisturbed flow velocity at the cross-section of the hole. The 
problem is thus essentially one of flow at shallow angles through a perforated 
wall, which can be likened to the flow about a series of wings at small angle of 
attack, the solid portions between the holes representing irregularly-shaped 
wings. The pressure difference set up by the perforated wall is then proportional 


: 
: 
| 
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to the lift on the solid portions and thus to the mean normal component of 
velocity, so that one is again led to the boundary condition (3) for perforated 
walls. 

This linear behavior is illustrated in Figures 1 and 2. Figure 1 shows 
experimental results for flow past a transverse slot and a circular hole, as well 
as the theoretical linear relationship of equation (3) and quadratic relationship 


Theoretical curves 


Linear law —s 


Quadratic law —--- 


10 


at es Ee 


Experimental points 


Circular hole fe) 


Transverse sloF A 


Figure 1 


Characteristics of a one inch transverse slot and a four square inch circular hole. 


of equation (4). Note that positive values of v,,, the mean normal velocity 
across the opening, represent an outflow through the wall (that is, the main 
flow is considered to be in the lower half plane, y < 0). The theoretical relation- 
ships are derived in more detail in [37]. Figure 2 shows experimental results for 
flow past perforated walls of different ratios 6 of open to total surface. Again 
positive values of v’ represent an outflow through the wall. While these experi- 
mental curves verify the linear relationship, it is seen that they actually pass 
slightly to the right of the origin, indicating that there is some outflow even 
where there is no pressure difference across the wall. This is caused by turbulent 
mixing through the wall, as explained in [34]. 

It should be pointed out that, in order to generate lift, a Kutta condition 
must be imposed on the flow at the trailing edges of the solid portions (leading 
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8=196% O 
6=22% A 
b=46% V 


Figure 2 


Wall characteristics of several types of perforated sheets. 


Figure3 
Perturbation velocity distribution along the wall for a single transverse slot. 
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edges of the holes). Thus the flow possesses a distinctive orientation, which 
leads to an unsymmetrical flow pattern even though the basic disturbance is 
symmetrical. Experimental results which bear out this contention are shown 
in Figure 3, where the theoretical flow past a simple transverse slot with a 
Kutta condition at the upstream edge, as calculated in [37], is compared with 
experiment for a value of 4P/q = 0-256. (The quantity g= 1/2 @ U? denotes 
the undisturbed dynamic pressure.) 


B. Mean Boundary Condition 


It is clear that the flow at some distance from a slotted wall which is large 
relative to the slot spacing may be considered as a uniform mean flow with a 
normal velocity component v’. For slots which are small and closely spaced 
this flow may then be identified near the wall with the flow in a slotted wall 
test section containing a model, thus yielding a mean boundary condition. 
Depending on the orientation of the main flow direction in the wind tunnel with 
respect to the slots, the cases of longitudinal, transverse, and swept slots will 
be obtained. Since, as indicated by Figures 1 and 2, the behavior of transverse 
slots and perforated walls is similar, the mean boundary condition for a perfo- 
rated wall may also be found. The treatment of these cases collectively rather 
than individually results in a much better insight into the behavior of partially 
open walls. 

The solution for a uniform incompressible flow with a mean velocity com- 
ponent v’ normal to an infinite number of slots of width 2 @ and spacing 2/ 
(Figure 4), over which the perturbation potential shall vanish, can be constructed 
as described in [37]. The complex conjugate perturbation velocity c*’ must pos- 
sess poles at the locations + a— 21, where 1 = —oo,..., — LO; 2 1s eee 
+ oo. The location of the zeros (stagnation points) shall be left unspecified and 


Pressure = 


beet at, 


Pressure=P 


Figure 4 
Flow past a slotted wall. 
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denoted by s — 2/. Then, in terms of the complex variable z= x +12: 


(ee s— 2nl) 
ised = A oo eta a = — 
|| x —a— 20) @+a—20)) (5) 


sally sin (2/2 1) (@—S) t 
/sin ( m/2 i) (z — a) sin (xf2 1) (2 + a) 


where A is a real constant. 
Since the real part of this function vanishes for y = 0 and 


2nl—@e =< 2nb+ 4, 


while the imaginary part is zero everywhere else on the x-axis, this solution 
satisfies the required boundary conditions. As long as s is not determined by 
some additional condition (Kutta condition or, in the case of longitudinal slots, 
requirement of symmetry) this represents an infinity of solutions. It contains 
all possible flows satisfying the stated boundary conditions for which the phy- 
sical half plane is mapped into a single Riemann surface of the hodograph. 
The two basic flows, a linear combination of which forms all these possible 
solutions, are obtained for s = 0 and s = J. . 
If the conjugate perturbation velocity is integrated and the integration © 
constant is determined such that the perturbation potential is zero within the 
slots and the perturbation stream function is zero on the solid portion of the 
wall, the complex perturbation potential is obtained as 


ga ey es [sin ws sin (x 2/2 1) +]/ sin®( 2/21) _ 1| 


sin (% o/2) sin? (at 0/2) 


(6) 


+ 4c0s a, in| cos (x 2/2 2) a |/ cos? (x 2/2 1) 1}} 


cos (x o/2) cos? ( 0/2) 


where o = a/l denotes the ratio of open to total wall area. 

The complex perturbation potential S’ and the conjugate perturbation 
velocity c*’ are now evaluated at some distance y = — Ay from the wall such 
that Ay/l > 1, i. e., where the disturbances caused by the wall are averaged. 
Note that the undisturbed flow is again located in the lower half plane. Clearly 
for sufficiently closely spaced slots this condition may be satisfied quite close 
to the wall, yielding the mean conditions 


iy = A {z ew is/21 Bik 7S 
+21 sin 21 eS a Insin 2 oo fe 
eee mS to 
+ 4 —— cos-5- Incos : a 
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es pa een 2 
ee ns aay. eee -ins/21 \ 
Cc dz U 1) = Av , (8) 
ino Wile Tiles ws v: . v's 
v =w tan a7 =Asin—=—.. (9) 


Separating real and imaginary parts in equation (7) and evaluating the 
constant A from equation (9) yields the perturbation potential at y = —Ay: 


Py= dy = = % Cot - ; —v Ay+v’ = Insin™* 
If this mean flow were extended to the wall at y = 0, it would have appro- 
ximately the potential ~,_9 = ~,= 4, +0’ Ay. Thus the mean boundary con- 


dition at the wall for an incompressible flow becomes 


Sekt uO 
=o x cot 5 Se Ce “<= 
e = En ee Insin a 


This result may be applied to slots swept back at an angle 2 by replacing 
x with x cos and may be further modified for use in linearized compressible 
flows with the aid of the Prandtl-Glauert rule. The eas then becomes 


@ V1 — M2, cos? Q =0' x cosQ cot + y — Se tata 2, 

This boundary condition is now matched with the actual flow in a slotted 
wall test section containing a model; that is g’ and v’ are identified with the 
perturbation quantities y’ and v’ induced by the model, and 0’ x is identified 
with the stream function — v, indicating the quantity of fluid leaving the 


test section. Then 


y’ Vi NE cost Le ap! cot (—F 7) cos 2 + v’ — Insin=” (10) 


21 2 


This is the mean boundary condition to be imposed at the upper wall for a 
linearized compressible flow inside a test section with a large number of slots 


swept back at an angle 2. 
In order to obtain a better physical understanding of the various terms, 


equation (10) is differentiated with respect to x. 


cos oY Lv‘) Insin =" (11) 


u Vi cos?,.2 = v’ cot 3 


If longitudinal slots are employed (cos 2 = 0) the first term on the nght 


hand side vanishes. Then 
“ue +Cv,=0 (12a) 


where 
C= Ine (12b) 
1 Z 
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The boundary condition for longitudinal slots may be interpreted as stating 
that the pressure difference across the wall is proportional to the streamline 
curvature and is not influenced by compressibility. Since the basic theory is 
identical to slender wing theory, the boundary condition will remain the 
same for slightly perturbed subsonic and, excluding waves of finite strength, 
transonic and supersonic flows. 

For the case where cos 2 + 0 a Kutta condition must be imposed on the 
flow at the upstream edge of the slot. Hence the perturbation velocity possesses _ 
a zero there, and the quantity s/ becomes —o, thereby determining the magni- 
tude of the coefficient of the first term on the right of equation (11). The 
second term on the right, containing as a factor the change of the normal velo- 
city over a single slot spacing, is then much smaller than the first term, since 
the ratio of the coefficients of the two terms is of order unity. Thus the second 
term can be neglected as soon as the Kutta condition takes effect. The boundary 
condition for a test section with non-longitudinal slots then becomes 


u'+ coté, v' =0 ~ (13ah 
where 


uO if 
cot, = cot; act ae (13 b) 
This boundary condition, which contains a Mach number effect, is valid for 
supersonic as well as subsonic flows, provided the slots are swept well behind | 
the Mach cone. 

A similar condition, with a different constant cot 8, however, will be correct 
for perforated walls. This constant will be smaller than for the case of transverse 
slots (2 = 0), and the Mach number effect will be greatly reduced if the per- 
forations are circular. 

It should be noted that the ratio of open to total wall area can be reduced 
in the swept slot case by increasing the angle of sweep and in the longitudinal 
slot case by decreasing the distance between slots (relative to test section 
height). This is of considerable importance, since the friction losses increase 
rapidly as the percentage of open wall area is increased. 


4. Subsonic Interference in Partially Open Test Sections 


A. Wall with Longitudinal Slots 


As already mentioned, the early work on subsonic interference in the slotted 
wall tunnel was carried out by using the exact boundary conditions that u’ = 0 
at the openings and v’ = 0 at the solid portions of the wall; it soon became 
apparent, however, that the use of a mean boundary condition furnished 
results of sufficient accuracy. In fact, since the work with the exact theory 
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must be carried out in some instances numerically with infinite series, of which 
only a finite number of terms can be considered, the use of the mean boundary 
condition, besides simplifying the problem considerably, seems to yield more 
accurate results. 

Subsonic interference, using the linearized approach, has been determined 
for a number of cases by various investigators using several different techniques. 
The most promising of these methods seems to involve the use of Fourier 
integrals in the two-dimensional case and integrals of Bessel functions in the 
three-dimensional case of a circular tunnel. In such a way blockage, buoyancy, 
angle of incidence, and streamline curvature interferences have been evaluated 
for tunnels containing a two- or three-dimensional doublet, a lifting wing of 
finite thickness, or a source simulating the wake of a model. (See particularly 
[20], [21], [22].) 

If a rectangular wind tunnel with solid side walls and slotted top and bottom 
walls is to be employed, the interference for two-dimensional models can be 
modified to obtain approximate values of interference for three-dimensional 
models in the tunnel. 

A rectangular wind tunnel with all four walls slotted will yield approxi- 
mately the same interference as a circular tunnel of equal cross-sectional area. 

The method to be followed in a tunnel with solid side walls and transonic 
(partially open) top and bottom walls containing a two-dimensional model will 
be discussed in some detail in a later section. 


B. Perforated and Porous Walls 


Since the mathematical formulation of the mean boundary condition at a 
perforated wall is the same as for a porous wall, the work done on the porous 
wall tunnel can be applied for perforated walls also. Again the method of 
Fourier integrals or integrals of Bessel functions will lead to the desired results 
in a most elegant fashion. Other techniques, such as Fourier series and image 
methods, have also been applied to the problem. 

The cases evaluated include the circular tunnel containing a three-dimen- 
sional doublet or lifting wing of negligible span and the tunnel with two per- 
forated and two solid walls containing two-dimensional singularities. (See [27], 
[28], [30].) The treatment of the finite wing in the latter tunnel encounters 
difficulties which have not been overcome as yet, since for the flow in the 
Trefftz plane (where w’ = 0), the boundary condition of equation (13a) becomes 
that of the closed tunnel (v’ = 0). Thus the calculations would have to be 
carried out in the complete three-dimensional system. Images of the lifting 
wing itself will be similar to the images of a two-dimensional vortex between 
perforated walls, while the images of the trailing vortices will be the same as 
if a solid boundary were employed. The angle of incidence interference due to 
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the trailing vortices will therefore be the same in a perforated test section as 
in a closed tunnel, since the presence of the solid lateral portions of the per- 
forated wall insures the existence of a Kutta condition and the resulting 
boundary condition of equation (13a). It may be possible to obtain the angle 
of incidence correction required in the tunnel with two perforated and two 
solid walls on this basis. 

Interference can be estimated for the rectangular tunnel with four per- 
forated walls by comparing it with that of the circular tunnel of the same 
cross-sectional area. 


 C. Interference in the Tunnel with Two Solid and Two Transonic Walls 


It has been shown that the exact boundary conditions for a wall with 
longitudinal slots and a wall with perforations or non-longitudinal slots may 
be approximated by the mean boundary conditions of equations (12a) and 
(13a) respectively. Thus the general mean boundary condition for a test section 
of height / with transonic top and bottom walls, assuming incompressible flow, 
may be expressed by y = +A/2: 


gy’ =+y’ cot#) + Cv’ (14) 


where the constants 9, and C are given by the following relations. 
Longitudinal slots: 


cotd; = 0.4 C= eed Insine= =; (15a) 
I 2 
transverse or swept slots: 
cote, = cote Cost) astral: (15b) 
perforated walls: 


cotO, = const > 0, C=0. (15c) 


The flow in a wind tunnel containing such walls together with solid side 
walls must satisfy this boundary condition at the partly open walls. If the 
tunnel contains a two-dimensional model the problem is entirely two-dimen- 
sional, and the complex perturbation potential of the flow can be expressed. 
as a Fourier integral. That is, 


®'(s) = Bi) + ©,,(2) = [A ot” a (16) 
where = 
A).= BO IG). (17) 
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In these relations ®;(z), B(t) and @,,,(z), I(t) refer to the complex potential and 
its Fourier transform for the basic free flight perturbation flow and the inter- 
ference flow (due to the influence of the walls) respectively. A(t), B(t) and I(2) 
are complex functions of the real variable ¢. 

Since the interference at the body is due to the effects of the walls as 
influenced by the body, the basic flow may be simply approximated (doublet 
for thickness of body, single vortex for lift, source for wake, etc.) for calculation 
of interference effects. The basic flow, then, is known and becomes 


co fo.) 
i (2) = a Bit) e'**!* dt = [{ BY) et 4 Bla) e-*"} at (18) 
—co 6 

Since the model is located on the x-axis at the center of the test section, 
the flow will possess singularities only there. Consequently the integrals of 
equation (18) must converge elsewhere, and it becomes necessary to distinguish 
between B(t) and B(—t) for positive and negative values of y respectively in 
order to insure such convergence. 


y>O, z= iph: BQ=L-{G, Gp}, B-H=0, | re 
y<0, z=-iph: BY) =0, B(-)=L-{G(-ip)}. | 


In these expressions # is some complex number whose real part is positive, 
and Lt g(p) }= 7(t) denotes the inverse Laplace transform of the function 
g(p) such that: 

g(p) = | L-*{g(6) }e* dt = L {1}. (20) 


0 


These inverse Laplace transforms can be obtained from tables for all flows of 
interest. 

The general mean boundary condition for partly open walls and an incom- 
pressible flow [equation (14)] may now be rewritten in terms of the complex 
perturbation potential: 


; . ’ dD’ dD’ \*) _ 
y= 28 et t% p’ is et t% @*' + C sin}, | -P ( a | = 0) (21) 
since: 
, - pe ao (d®’ \* 
29’ =@' + OD’, 2iy' = @' — o*’, Oe az - (=) 


' where the asterisk denotes the conjugate of the function. If equation (16) is 
substituted into this general boundary condition and use is made of the fact 
that the resulting equation must be satisfied at all points on the boundary 
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(i. e., for all x), the condition to be satisfied by the coefficient function A(t) is 


found. 
For y = h/2 (top wall): 


A(t) a* (—t) = A*(—?2) ad) , (22a) 
for y = — h/2 (bottom wall): 
A(t) a(t) = A*(—2) a*(—2) , (22b) 


where: 


a(é) = el? (at +4 me sino) f (23) 


If the result of equation (17) together with the values for B(t) and B(—?) 
from equation (19) are now substituted into these equations and the resulting 
system is solved for J(¢) and J(—?) one obtains: 


LO = sag erty LEE BH) a(t) + LB (i Jal], (24 


ae ; 

I( d= Sag y LEO to, (¢ ph) hax (t) )+ L- 14 Dy ( (—i ph) }a(—t - (25) | 
Thus the eee interference potential may be found in terms of the 

integral of a function involving known inverse Laplace transforms, since 


= (9 ei tlh ae i t) Caney di (26) 


Generally the integration can be carried out in closed form only if C = 0, 
i.e., if perforated walls or non-longitudinal slots are employed. However, it is not 
too difficult to perform the integration numerically for the derivative, d®,,,/dz 
at z = 0, for C + 0, thus obtaining the interference velocities at the origin. 

As an example, the complex interference potential due to a two-dimensional 
doublet located at the center of the test section shall be evaluated. The inverse 
Laplace transforms for the free flight complex perturbation potential of the 
doublet O; = M/2 zz are: 


y>0: LAO Gph)}= i“ 1H), (27a) 
y<0: L{G,(-iph)}=i 5 1 (27b) 


where 1(¢) denotes the Heaviside unit function. 
1) =0 <0), (28a) 
1=1 6S On (28b) 
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By combining the results of equations (24)-(28b) the complex interference 


potential is obtained: 


®,,(2) =i 5 lee Ae tt ee 1] eth ae (20) 


For longitudinal slots 6, = 2/2, and hence 


a(t) =sel@ (14 © t), (30) 


ees aes RIC N39 a) “a 
0 
which at z = 0 yields the conjugate interference velocity as a function of C/h: 


(Ct),-0= (“), = sp /t toa ae 


1 — (Cth) 


Since C/h varies from zero for the completely open test section to infinity 
for the completely closed test section, it is more convenient to plot the inter- 
ference velocity as a function of the parameter 
1 — (C/h) 


bce bh eG: 


which varies from —1 for the closed test section to +1 for the open test section. 
This tunnel parameter TJ is plotted in Figure 5 in terms of o, the ratio of open 


a ee Sanh 
enews fe ser Ae | 
Bae he Be ee 
(= SE ee) EA ee 


Figure 5 
Tunnel parameter T for various slotted wall test sections. 
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to total wall area, and h/2/, the ratio of tunnel height to slot spacing. Within. 


linearized theory T is not influenced by compressibility. 

Since the integral of equation (32) is real, the angle of incidence interference 
v;, Vanishes at z = 0, and hence the relation directly yields the blockage inter- 
ference u,,, which is plotted in Figure 6 as a function of T. It is seen that for 
T = 0:25 (C/h = 0-60) the blockage interference vanishes. 


Figure 6 
Blockage interference onatwo-dimensional doublet located at the center of a slotted wall test section. 


For non-longitudinal slots or perforated walls C = 0, so 
a(t) = elt/2) +79, 


For this case, with the aid of [44], the integration of equation (29) can be 
carried out, yielding an interference potential: 


Nhe M 20) 2/h TZ M 
®,,= 57-6 {coth eae 1} — (33) 
Since this function is analytic at z= 0, it may be expanded'i in a Taylor 


series about this point, thus yielding the various interferences at “the location 
of the model. 


a Pe 
Din Fy ee ete 


(34) 


463 — 6m 02 + 2°22 0, ( 2\2 
= ae (=) +o}. 


The blockage and buoyancy interferences at z = 0 are then obtained as: 


(Winano= A( GE"), = saz (2 68-2054 4), (35) 


~~ 
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OUin = a’®@,, M / 4 o 2 = 
(i) Ri di ine 7 Pots, 9% —- 29+ 3 zt) (36) 


where “ denotes the real part of the function. 

Making use of equation (35), and-recalling that cot 6) is a positive constant 
[equations (15b) and (15c)], and hence 0) < 2/2, it is seen that for a wind 
tunnel with perforated walls or non-longitudinal slots the blockage interference 
for a two-dimensional doublet located at the center of the test section vanishes 


at the doublet for 
Es 1 
6.= ar (1 _ a (37) 


However, if the wall is adjusted to meet this condition, a buoyancy correc- 
tion, given by 


, 


Ou; n M Gi 
Coorg Mere poe oye Ge) 
must still be applied. This is one of the distinguishing features between the 
perforated wall and longitudinal slots,.where, for the doublet, the second 
derivative of the complex interference potential, and hence the buoyancy cor- 
rection, vanishes at z = 0, as is easily shown from equation (31). 

If the Prandtl-Glauert rule is applied to the flow in the perforated wall 
tunnel, it is seen that the required boundary condition is Mach number de- 
pendent. That is, 


| i oS 1 tan cot 0, 39 


where the subscripts c and 7 refer to the compressible and incompressible cases 
respectively. 

A comparison of equation (39), which is also valid for a wall with non- 
longitudinal slots, with the boundary condition for compressible flow past the 
upper wall, as expressed by equations (13a) and (13b), shows that 


Ce | ee ae 
cot@, =icot a we — Me, * 


Combining this with the condition for vanishing blockage interference, as 
expressed by equation (37), yields 


= 2 coe {AME tan -7-| “ 
girs Col i= mM, 1 ys (40) 


as the condition for zero blockage interference on a doublet at the center of a 
test section with non-longitudinal slots for compressible flow. For transverse 
slots (Q = 0) the Mach number influence vanishes and 6 = 42:3% is obtained. 
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D. Discussion of Some Additional Interference Results 


Making use of the method just described a variety of results have been 
obtained for wind tunnels with two transonic and two solid walls. (See) for 
example, [20].) Figure 7 shows, for incompressible flow, the variation of the 
blockage interference along the axis of a wind tunnel for longitudinal slots and 
perforated walls due to a doublet located at the center of the test section. The 
non-symmetrical behavior for the perforated wall tunnel, which is due to the 
Kutta condition employed at the wall, is clearly evident. 
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Figure 7 


Blockage interference along the tunnel axis due to a two-dimensional doublet located at the center 
of a partially open test section. 


Figure 8 illustrates, for incompressible flow, the variation of blockage inter- 
ference along the axis in a tunnel with longitudinal slots due to a two-dimen- 
sional source of strength Q which is located at the center of the test section. 


This interference vanishes at the source for all cases except the completely | 


closed tunnel. It can be shown [20] that the buoyancy correction due to the 
source vanishes for the same value of T for which the two-dimensional doublet 
yields zero blockage interference. Thus in the tunnel with longitudinal slots 


blockage and buoyancy interference due to both the model and its wake vanish 
if the correct wall is chosen. 


By considering the flow in the Trefftz plane the angle of incidence correction | 


can be determined for longitudinally slotted walls. For the two-dimensional 
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wing in such a tunnel this correction obeys the extremely simple relation [20]. 


Ha 
U 4h 1+ 2 (Ch) 


(41) 


where d denotes the wing chord, C, the lift coefficient, and C the tunnel para- 
meter as defined by equation (15a). 


50 


ee] pe 
on 


- 


ca ee 


li 
i 


Oo 
+ 
oO 
& 
>|x 


Figure 8 


Blockage interference along the tunnel axis due to a two-dimensional source located at the center 
of a slotted wall test section. 


The angle of incidence interference for incompressible flow past a finite wing 
at the center of a square test section with longitudinally slotted top and bottom 
walls, as derived in [20], is given in Figure 9 for different values of 8, the ratio 
of effective wing span to tunnel width. The cross-section area of the test section 
is denoted by A and the wing planform area by S. Figure 10 gives the value of 
this correction for test sections with different values of A, the ratio of tunnel 
width to height, for wings of small . 

Using an approximate method described in [20], the blockage interference 
at z = 0 for incompressible flow past a finite wing at the center of a square 
test section with longitudinally slotted upper and lower walls may be deter- 
mined. These results are plotted in Figure 11, where M denotes the strength 
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Figure 9 


Angle of incidence interference on a finite wing located at the center of a square slotted wall test 
section. 


Figure 10 


Angle of incidence interference on a finite wing of small span located at the center of a slotted wall 
test section. 
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Figure 11 
Blockage interference on a finite wing located at the center of a square slotted wall test section. 
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Figure 12 


T as a function of f for zero interference on a finite wing located at the center of a square slotted 
wall test section. 
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of the three-dimensional doublet distribution which represents the wing. Note. 
that the value of 8 = 0 would correspond to a simple three-dimensional doublet 
(body of revolution). 

From Figures 9 and 11 it can be seen that for each value of f there is a 
particular value of T for which the interference vanishes. These values are 
plotted in Figure 12 for a square test section. It can be seen that it is impossible 
to find a wall configuration for which both interferences vanish. 

Figures 13 and 14 give experimentally determined values of the pressure 
coefficient at two points on the surfaces of two geometrically similar half 
diamond airfoils. The airfoil is located near the lower solid wall of a 9” x 9” 
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Figure 13 
Pressure coefficient as a function of Mach number at pressure tap 1 on a half diamond airfoil. A 
sketch of the model is given in the insert. The chord lengths d are 5-5 and 3-5 in. for airfoils a and b 
respectively. 
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test section, as indicated in Figure 13, and the results were determined by 
using both solid and slotted upper walls. The exact experimental setup is de- 
scribed in more detail in [14]. While all of the experimental points have been 
plotted, to avoid the confusion of having many curves very close together, 
not all of the curves have actually beén drawn. The curves labelled ‘closed 
corrected’ were determined by modifying the results obtained in the closed 
test section by the use of a correction factor, as described in [14], and approxi- 
mately represent free flight conditions. The deblocking effect and the reduction 
of interference for the slotted wall test section are obvious. 


5. Supersonic Interference in Partially Open Test Sections 


By following the analysis of the previous section it is possible to determine 
the wall interference acting on a model in a wind tunnel, provided that the flow 
is subsonic. If, however, the flow is locally supersonic, shock waves will form 
which, in general, will be reflected from the tunnel walls. The interaction of 
these reflected waves with the model will introduce additional interference 
effects which must be calculated. This problem has been studied by various 
investigators for both slotted wall and perforated wall tunnels (see, for 
example, [38]-[43]). 

As was previously pointed out, the constant of proportionality contained 
in the mean boundary condition for longitudinal slots is not influenced by 
compressibility. Hence, if waves of finite strength are excluded, this boundary 
condition [as given by equations (14) and (15a)] is directly applicable to linea- - 
rized supersonic flows. The boundary condition for the perforated wall is, 
on the other hand, Mach number dependent, the constant of proportionality 
experiencing a rapid change with Mach number. 

Since the exact problem of determining the wall configuration for which the 
incident waves will not be reflected involves the study of wave-boundary layer 
interaction, only a few simple examples assuming potential flow throughout 
the flow field shall be considered here. This will be sufficient to obtain some 
insight into the nature of the problem. 


A. A Weak Two-Dimensional Wave Incident on a Wall With Longitudinal Slots 


Consider the problem of supersonic flow in a test section containing a two- 
dimensional slender body. The inclined portion of the body will, within linea- 
rized theory, generate parallel isentropic waves which will be reflected by the 
tunnel walls. These reflected waves may reach the body, causing errors in the 
pressure distribution. 

More specifically, consider a symmetrical body of length / and thickness 
distribution t(x) = 2 y9(x) at zero angle of attack located at the center of the 
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S Ai 


0-5 0-6 07 
Figure 14 


Pressure coefficient as a function of Mach number at pressure tap 2 on a half diamond airfoil. 
Symbols are the same as those used in Figure 13. 


test section (y = — h/2), as shown in Figure 15. By reason of symmetry, then, 
it will only be necessary to consider the flow over the upper half of the body. 
For convenience, use will be made of the notation: 


S=%— ViCOtK, 1 = 4 cote 


where « = sin-11/M,, denotes the Mach angle. The perturbation potential 
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Upper wal! 


Acoke 


Z 
Figure 15 


Supersonic flow through a test section containing a two-dimensional slender body, showing the 
associated wave pattern. 


for such a flow will then be given by the relation 


p = f(s) + feld) (42) 


where the function /,, which is constant along the line s = const, represents 
the incident disturbance generated by the body, and /,, which is constant 
along the line ¢ = const, represents the disturbance reflected from the wall. 
Since the body is located in the interval 0 < s < J, this perturbation potential 
will vanish for s < 0, thus insuring that /,(0) = /,(0) = 0. 

The boundary condition at the upper wall is given by 


wa) Bata x Cm 0. (43) 


If the potential of equation (42) is now inserted into this boundary condition 
and use is made of the fact that /,(0) = f,(0) = 0, the function /,(x),_9 may 
be obtained in terms of /,(),_9; that is, the functions are related at the wall. 
By moving along the two lines s = const and ¢= const passing through the 
common point s = ¢ = x on the wall, it is seen that 


fi(*)y-0 = h(s) =fh@), fel*)y-0 = held) 


and thus the solution for /,(¢) becomes 
t 


fal) 4 AG) 20 ne greek | AU) ettanalC gy (44) 


0 
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If the integral in this equation is integrated by parts and the resulting 
equation is then differentiated, a relation between the derivatives of the 
functions /, and f, is found. 

t 


df(t) _ a(t) _ 7 tan | ttana/c i ee gttanal/C gy (45) 


Sat dt cS 
0 


By differentiating the potential of equation (42) with respect to y, the 
vertical component of the velocity is found as 


of a OP meaty ae pu j. (46) 


The boundary condition that the flow be tangent to the body may be 
expressed as ts 
ee & ind 
eh ues = | (47) 
Since, on the body d/,(¢)/dt = 0 in the interval —h cota/2<%<h cota/2, 
a combination of equations (46) and (47) will yield, for this region 


af, (s) — —_U tang (o>) 
ds ax }*=s-—(heota/2) 


(48) 


The value of dj,(t)/dt may now be found on the body in the interval 
hcota/2<x<3hcota/2 (i.e., due to the first reflection) by substituting the 
result of equation (48) into equation (45). 


eS Ut ( ax beter ) 


42U se gnttnale [ (Ys) ef tana/C dt. 
X | x=t—(hcot a/2) 
0 


The horizontal component of the perturbation velocity in this interval is 
given by 


Oe a eis) a ee aa U tana (oe 


U 
Ox ds of dt ax | ase dol 


(50) 


where the second result is obtained with the aid of equations (46) and (47). This 
result, together with the value of d/,(t)/dt from equation (49), is valid for the 
body in the test section. The value of u’ for the body in free flight is given by 
letting d/9(¢)/dt = 0 in equation (50). The blockage interference, which is the 


difference between these two values, then becomes . 
5 4 
7 cota <i e Ricoto. ae =a ara? (51) 


x 


: 
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It is seen, then, that the interference on the hody due to the first reflection 
is twice the value of the disturbance reflected from the wall, as given by 
equation (49), 

The interference may be determined in a similar manner for x > 3 h cota/2 
(i. e. , for successive reflections of the disturbance) by use of equation (45), 
with the exception that the value of d/,(¢)/d¢ will no longer be given by equation 
(48) but must be determined from equation (46) with the aid of the previously 
calculated values of d/,(t) /dt. 

As an example of the application of this result, consider the flow past the 
body of Figure 16. For such a body 


E aa 

dVq ae 
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Since the slope of the body is zero for x << —hcota/2, no disturbance is 
created in this region, and hence the blockage interference will be zero for 
x < hcota/2 (¢< 0). A combination of equations (49) and (51) under the 
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Blockage interference due to the reflection of a disturbance from a slotted wall. A sketch of the 
model and the associated wave pattern is given in the insert. 
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conditions of equations (52) then gives the interference for 


cota<*< - hceota (O<t<hcota) 


as 
| 2.U 6 tana (Lim 2 ee ea 0 eee | 


= (53) 
| 4U dé tang etna (Avenel 4)" (te <a cota) | 


U; n(t) 


These results are plotted in Figure 16 for different values of the parameter 
Ax tana/C. The incident wave, which is initially reflected as if a solid wall 
were present, is rapidly followed by waves of the opposite sign, the effect of 
which is to reduce the interference on the body as indicated. 

It is possible that in the actual physical case, where the waves are not 
parallel, but converge slightly, some wave cancellation may occur before the 
reflected waves are able to disturb the flow on the body. However, an investi- 
gation into wave-boundary layer interaction must precede any definite state- 
ment in this respect. While in one case excellent cancellation may be obtained. 
by removing the boundary layer through a converging wall with a large number 
of slots, conditions may be very different for the case where only a small 
number of slots are used and the walls are diverged sufficiently so that the 
boundary layer is not removed. 


B. A Weak Two-Dimensional Wave Incident on a Perforated Wall 


The boundary condition at the perforated upper wall of the test section, 
as given by equation (13a), may be rewritten as 


y=0: y, + coté, vy, = 0 (54) 


where the constant of proportionality, cot 6,, is Mach number dependent. 
If the supersonic perturbation potential of equation (42) is now inserted into- 
this boundary condition, the condition at the wall becomes 


af, (4%) _ coté, — tang [ df,(#) 
| dx ihe ce et ce laa See 


It is seen that a definite attenuation of the disturbance takes place. In 


particular, the incident disturbance is completely cancelled by the wall for the: 


condition 4 


cot = tang = =e 
(Mz — 1 


(56): 
Thus, by proper choice of the constant 6,, the boundary condition for per- 


forated walls can be made to simulate the free flight streamline exactly for a 
particular Mach number, within linearized two-dimensional theory. 
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C. Three-Dimensional Waves 


If longitudinal slots are employed ina test section containing a three-dimen- 
sional model, the waves reflected from. different portions of the wall may 
interact and thereby be attenuated considerably. Thus the two-dimensional 
analysis already carried out yields results which are probably only a very 
conservative estimate of the actual three-dimensional performance. 

For a three-dimensional model in a perforated wall test section, however, 
the performance will not be as good as for the two-dimensional case already 
discussed, since the condition given by equation (55), which made possible the 
elimination of blockage interference-by proper choice of the wall configuration, 
is no longer satisfied. 

Thus, if a comparison is made of test sections with perforated walls and 
walls containing longitudinal slots, where the grain sizes (slot widths and hole 
sizes) of the two boundaries are comparable, and the subsonic portion of the 
boundary layer is removed in both cases, the advantages of perforated walls 
for two-dimensional flows at high transonic speeds are obvious. For three- 
dimensional flows, however, it seems that the perforated wall becomes less 
effective and the slotted wall more effective in reducing interference effects. 


6. Conclusions 


It has been shown that for subsonic flow in a test section with longitudinal 
slots the blockage interference on the model can be eliminated by proper choice 
of the slot configuration. It has also been pointed out that the buoyancy inter- 
ference on the model is eliminated by the presence of the slots, regardless of 
slot size. In [20] it is shown that the blockage interference on the wake is also 
eliminated by the presence slots of any size, and the buoyancy interference on 
the wake will be eliminated for the same slot configuration as is required for 
the elimination of blockage interference on the model. Thus these four inter- 
ferences can be completely eliminated by use of slotted walls. The angle of 
incidence correction has been evaluated in [20] for a rectangular test section 
with top and bottom walls slotted and in [21] for a slotted circular test section. 
The streamline curvature interference is evaluated in [21] for the slotted circular 
test section and in [22] for the rectangular test section with top and bottom 
walls slotted. 

It has also been shown that it is possible to choose a rectangular test section 
with top and bottom walls perforated which, for subsonic flow, will eliminate 
the blockage interference due to the model. Such a test section will also insure 
the elimination of the buoyancy interference due to the wake. The buoyancy 
interference due to the model and the blockage interference due to the wake, 
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contrary to the slotted wall case, are not eliminated by the perforated walls. 
This presents no problem, however, since the buoyancy correction for the model 
is determined by equation (38), and the blockage correction for the wake may 
be similarly easily determined. At the present time the angle of incidence inter- 
ference has been evaluated only for the circular perforated test section [28] and 
the streamline curvature interference has been evaluated for neither the rec- 
tangular nor the circular test section. Thus, while the interference problem has 
been solved for the slotted wall test section, much work still remains to be done 
for perforated walls. 

For a two-dimensional supersonic flow through a test section with a large 
number of longitudinal slots there is no initial attenuation of the disturbance; 
however, this initial undiminished reflected wave is followed by reflected 
waves of opposite sign, thus reducing the interference on the model. Due to 
wave interaction, better attenuation is expected for the three-dimensional 
tunnel than for the two-dimensional tunnel, except for the case of the circular 
tunnel containing a body of revolution, which exhibits a two-dimensional 
behavior. 

Contrary to the behavior of the slotted wall, the perforated wall seems to 
perform more satisfactorily in the two-dimensional case than in the three- 
dimensional case. Whereas complete wave cancellation is possible for the two- 
dimensional-problem, such a result cannot be achieved for three-dimensional 
flows, since equation (55) is no longer valid. 

These interference effects are summarized in Tables 1 and 2 (opposite page) 
for convenient reference. 

Aside from the elimination or determination of interference for partially 
open walls, which has been discussed at length, there are other important 
considerations with respect to such walls which should be mentioned. One of 
these considerations is the establishment of supersonic speeds in the test 
section while using a sonic nozzle. This may be done either by confining the 
wall openings to the test section itself and using auxiliary suction or by exten- 
ding the openings into the diffuser walls, in which case no auxiliary suction 
equipment is necessary. Both of these methods are, of course, practical only 
for low supersonic speeds, and both require considerable power. In the former 
case the diffuser losses may be reduced by increasing the suction, so that some 
optimum operating condition may be found. 

In general it seems that a slotted wall test section will require more power 
than will one with perforated walls of the same grain size. For both of these 
cases the configuration of the wall openings at the leading edge of the test 
section must be modified in order to prevent an overexpansion there. . 

While these considerations have been only briefly mentioned here, they 


are, nevertheless, very important, and for a more detailed discussion the reader . 
is referred to [23]. : 
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Table 1 
Subsonic Flow 


Type of interference Walls with longitudinal slots Perforated walls 
Blockage, model May be eliminated May be eliminated 
Blockage, wake None May be evaluated 
Buoyancy, model None May be evaluated 
Buoyancy, wake May be eliminated May be eliminated 
Angle of incidence May be evaluated Known only for circular 

tunnel 
Streamline curvature May be evaluated Presently unknown 


Table 2 
Supersonic Flow 


Type of flow | Walls with longitudinal slots Perforated walls 


Two-dimensional No initial attenuation. Ini- | Complete wave cancellation 


tial reflection followed by | possible. 
waves of opposite sign, re- 
ducing the interference. 


Three-dimensional Better attenuation than for | Complete wave cancellation 
two-dimensional flow, ex- | impossible. 
cept for circular tunnel. 


(1] 


[2 


Lo. 


[3] 


[4] 
[5] 


[6] 
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Zusammenfassung 


Um von Windkanalversuchen auf die Verhaltnisse beim freien Flug schliessen — 


zu kénnen, miissen die Einfliisse der Kanalwande auf die Geschwindigkeitsver- 


teilung am Modell beriicksichtigt werden. Bei schallnaher Str6émung sind diese ~ 


Einfliisse an Modellen tiblicher Grésse so stark, dass sie die Messungen voll-— 


standig verfalschen. 
Die hier beschriebenen Methoden beruhen auf der Berechnung von 4 Arten 


von Interferenzen in einem Punkt des Modells. Von diesen werden auf einfache ~ 
Art die Interferenzgeschwindigkeiten in den anderen Punkten des Modells abge- 


leitet. 

Ausgehend von einer gemittelten Grenzbedingung werden die Interferenz- 
geschwindigkeiten fiir teilweise offene Kanalwande berechnet. Durch die Wahl 
verschiedener Wandkonstruktionen kénnen einzelne Interferenzen eliminiert 


werden. Die Tabellen 1 und 2 zeigen eine Zusammenstellung der Méglichkeiten © 


fiir verschiedene Arten von teilweise offenen Kanalwanden. 


(Received: September 8, 1955.) 
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On the Structure of Tce-Forming Nuclei’) 


By RaymunD SANGER, Ziirich?) 


The cause of icing nucleability has become a more important question, 
since SCHAEFER and VONNEGUT [1], [2]3) discovered that silver iodide (Agl) 
can be used as a seeding substance to form ice. At that time the only other 
substance found to be similarly active was lead iodide (PbI,). The fact that they 
both, AgI at least to a considerable extent, display a hexagonal crystal like 
that of ice, led to the idea that the reason for the high degree of icing nuclea- 
bility in AgI must lie in its crystalline structure. This supposition was strength- 
ened by the further fact that the lattice constants of the three substances in 
question are virtually the same. 

Subsequently, however, the same research group observed a considerable 
degree of icing nucleability in certain other substances, whose crystalline 
symmetry bears no relationship either to that of ice or that of the hexagonal 
B-AglI crystal. And not surprisingly the lattice constants in these new cases 
also vary considerably from those of ice. Hence it is clear that neither the 
macroscopic crystal symmetry nor the lattice constants can tell us much 
about the icing nucleability of a substance. This truth has been confirmed by 
more recent experiments [3]. It has proved possible to measure separately the 
icing nucleability of the two variants of AgI which exist at normal tempera- 
tures: the hexagonal f-crystal and the cubic y-crystal. The measurement was 
the same in each case. This discovery gives us a pointer for further inquiry into 
the question. For we should notice that in both variants of AgI the coordination 
groups of the crystals are identical, and constitute in fact tetrahedrons. In the 
corners of each we find alternatively silver and iodine atoms. Now, the coor- 
dination groups of ice also constitute tetrahedrons, with very nearly the same 
geometrical dimensions. The only difference is that in the corners we find only 
oxygen atoms, while the hydrogen atoms occupy an intermediary position. In 
the most important two-dimensional surface planes of both the cubic and the 
hexagonal Agl crystals the tetrahedron edges (running alternately right and 
left) are arranged in exactly the same way as in ice. This makes it easier to 
understand why the same degree of icing nucleability was observed in each 
variant; and it also suggests that the cause of nucleability should be sought in 


1) Paper presented at the Tucson Conference on The Scientific Basis for Weather Modification 
Studies, Institute of Atmospheric Physics, University of Arizona, April 10-12, 1956. 

2) Physikalisches Institut der ETH. 

3) Numbers in brackets refer to References, page 217. 
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the surface forces. That is to say: in the forces which operate between the 


surface atoms of the seeding substance and the surface atoms of the water — 


(whether solid or liquid). 
WEICKMANN [4] has shown that there should probably be close connection 


between the phenomena of epitaxy and icing nucleability. In his treatment of 


the subject he mentions RoyER and FRIEDEL [5], [6], who state that variations 
of up to 12% in the geometrical spacing of the atoms are admissible during the 
orientated overgrowth of guest crystals on the sliced surface of a host crystal. 
It should not be forgotten, however, that in the examples given by WEICKMANN 
(for instance, the growth of ice crystals on mica) he is concerned with epitaxy 


in its usual or ‘classical’ form ; that is to say, when the guest crystal is allowed an — 
unlimited period of time to complete the process. By contrast, where supercooled — 


clouds are seeded with suitable icing substances, freezing takes place extremely 


quickly. We must therefore be careful not to make any comparisons between ~ 


the two phenomena without allowing due consideration to the time factor. 
WEYL has tried to explain the icing nucleability of AgI in an entirely 
different fashion [7]. Basing his argument on the knowledge that the ionic 
crystal of sodium chloride (NaCl) has an electric double layer on its surface, 
he attributes a similar double layer to the AgI crystal, in which the chemical 
binding forces are strikingly homopolar. In the case of NaCl the explanation 


for this double layer clearly lies in the considerably greater polarisability of the — 


negative chloride anion (Cl-) over the positive sodium cation (Na+). The. 


equilibrium of the interatomic forces on the surface of the crystal can then only 


— 


be explained in one way: that the cations in the outer layers have been forced 
relative to the anions, towards the inside of the crystal. Now, despite the — 


homopolar character of AgI, the same marked difference exists between the 
polarisability of the anion-like iodine and the cation-like silver. Thus, according 
to WEYL, an AglI particle must also possess an electric double layer on its 


surface, of which the negative side lies outermost (as. with NaCl). And this is 


of importance for any explanation of icing nucleability. 

For water droplets likewise possess an electric double layer, with again the 
negative on the outside. The reason for this is that the electric dipoles of the 
water molecules situated near the surface are, on the whole, orientated in- 
wards; that is, each parallel to its immediate neighbours. Whereas inside the 
droplet they are so arbitrarily disposed in all directions that they neutralize one 
another. This view of the matter is compatible with the following experimental 
finding: that on a statistical basis the fundamental arrangement of oxygen 
atoms in liquid water is tetrahedral, as in ice. And that in fact two kinds 


of tetrahedral structure complement one another: a $-tridymite-like structure 


(similar to that of ice) and a quartz-like structure. Below the point of maximum 
density for water, at 4°C, the ice-type of structure predominates. When the 
temperature falls below zero, then (according to WEYL) the surface layer, 
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with its electric dipoles turned inwards, at first hinders freezing. Only when 
the temperature falls far enough below zero, or alternatively when the double 
layer is disturbed by some external influence, can the transition from liquid 
to solid be initiated. 

Such a disturbance takes place, for example, when an Agl particle ap- 
proaches a supercooled water droplet. For the electric moment of the double 
layer of the seeding particle is itself orientated inwards and therefore, at the 
point of approach, points in a direction exactly contrary to the double layer 
of the water drop. The result is then that the AgI destroys the alignment of the 
electric dipoles of the water molecules. It is true that WEYL assumes that the 
AglI particle must penetrate the surface of the water. But this would not 
actually be necessary to produce a loosening in the orientation of the water 
dipoles. Indeed, the extremely rapid completion of the freezing process, when 
supercooled clouds are seeded with AglI, suggests rather that the mere proximity 
of its double surface layer is sufficient. 

Before we attempt to correlate, or even to compare, WEYL’s explanation 
for the icing nucleability of AgI with the theory of epitaxy, we must define 
more closely the preconditions necessary for any regular overgrowth of guest 
crystals on a host crystal. In the case of ice, the following features favour the 
successful development of epitaxy: 

1. That the coordination groups are tetrahedral, regardless of whether they 
are built up into a hexagonal or cubic lattice; 

2. that there is a uniform type of chemical cohesion; that is, a homopolar 
link like that of ice, with a similarly shaped potential field set up by the mutual 
chemical forces; 

3. that there is the same arrangement of the basic tetrahedrons in the 
surface layer as we find in ice, and with similar geometrical dimensions. 

The longer the time allowed for epitaxy, the greater the discrepancies which 
may be admitted between the geometrical spacings of the atoms in the seeding 
substance and the atomic spacings in the ice. The limit of possible difference 
will depend upon the degree to which the chemical binding forces are suited to 
one another. The kind of limit proposed by Royer and FRIEDEL on a simple 
percentage basis must therefore be interpreted with great caution, as we have 
already pointed out; especially if it is intended to draw from it quantitative 
interferences as to the icing nucleability of the substances in question. 

If epitaxy is to be appealed to at all, for an explanation of the rapid 
freezing of supercooled clouds when seeded with special nucleable substances, 
then all of the above mentioned conditions would have to be satisfied to the 
fullest possible extent. Where this is so we can understand why Ag (6 and y), 
cupric iodide (Cul), mercury telluride (HgTe), calcium selenide (CaSe) and 
lead iodide (PbI,) are such effective nucleable substances*). But the reason why 


4) For more recent measurements and results, see [3], [8], [9], [10], [12]. 
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cupric sulphide (CuS) and cupric selenide (CuSe) are also nucleable is less easy 
to see. Of course, the crystalline structure of CuS does contain tetrahedral 


units in layer formation, andthe tetrahedrons are arranged in the same way as~ 


in ice [11]. But on the other hand, the interatomic intervals do not correspond 
at all well to those of ice. Further, the immediacy with which all these seeding 
substances initiate freezing leads one to suspect that the mere optimal fulfilment 
of the conditions necessary for the epitactial growth of ice on an ice-forming 
nucleus is in itself not sufficient to render a plausible explanation of the me- 


chanism of freezing. 
If all the binary compounds, which have so far been found to be nucleable, 


are tabulated on the basis of the Periodic System, then we see clearly that one © 
of the two elements falls constantly on the left, or positive, side of the table, 


while the other without exception falls on the right, or negative, side [12]. In 


r 


addition, the element on the right always possesses a greater polarisability — 


than the element on the left. When we further ascertain that in practice all 
these compounds are insoluble in water), we arrive at last at something like 
a complete picture of the substances that provide ice-forming nuclei. And it is 
perhaps worth noticing in passing that this picture is strikingly similar to the 
list of semi-conductors. : 


We find, then, that active seeding substances possess a noticeable degree — 


of polarity, despite the predominantly homopolar character of their chemical 


bond. The greater polarisability of the anion-like element over the cation-like 


element causes an electric double layer to be formed on the surface of the seed- 


- 


ing substance — in the same way as that posited by WeEyL for AgI, and with 


the same orientation of the negative towards the outside. Are we, then, in a 
position to assume that WEYL’s mechanism for the freezing process of super- 
cooled clouds is in fact typical for all nucleable seeding substances? 


The essential point of WEYL’s mechanism rests upon the disorientation of 


the electric dipoles on the surface of the supercooled water droplet. But is it 
possible to accept this as a full explanation of all aspects of the freezing process? 
The problem of epitaxy cannot be ignored. And the important question arises, 
whether some combination of both mechanisms might not be possible which 
would bring us closer to an understanding of what actually happens. On the 
basis of such a composite picture we should expect that: despite the optimal 
fulfilment of all the conditions necessary for epitaxy, the orientation of the 
molecular dipoles on the surface of the supercooled water drop would at first 
prevent the water molecules from settling on the surface of the seeding substance 
in conformity with the potential pattern set up by the chemical binding forces. 


5 
* ) The substances such as AgNOg, KI, Nal, which are listened in the work by PRUPPACHER and 
ANGER [3], [12] as soluble in water, were atomized in the freezer by a fume process, since they 


could not be disintegrated mechanically. The resul i i 
: t 
Pasa aileslianle: s obtained in such cases cannot therefore be 
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The electric double layer, however, produced by the difference in polarisability 
of the atoms, breaks up the orientation of the water dipoles. And as the water 
molecules are then no longer hindered from settling in due form on the surface 
of the seeding substance, or from fitting themselves into an ice lattice, the 
freezing process can then be initiated-very rapidly. 

The experimental knowledge at present available is, however, so limited 
that it does not allow us to judge with even a fair degree of accuracy whether 
these explanatory mechanisms are correct. We must considerably increase the 
extent of our material observations along the lines of our present understanding. 
Above all we need more precise information concerning the initial temperatures 
at which nuclei become iceforming; secondly, we need to investigate a much 
larger number of substances — beginning preferably with binary semiconductors. 
And thirdly, attempts should also be made to further our knowledge of the in- 
teraction between the seeding substance and the water drop. This might be done 
by disturbing experimentally the surface of the ice forming nuclei. 

Finally it ought to be pointed out that, in the mechanisms discussed here, 
the seeding particle has been considered solely in the light of a freezing nucleus. 
The question whether there may not exist some form of sublimation nuclei, 
has not been touched on at all. And in conclusion we should not overlook the 
possibility of still further mechanisms which would also be capable of initiating 
the process of freezing. 


It is my pleasant duty to thank Dr. W. Epprecut for his very interesting 
contributions to the discussion. 


REFERENCES 


[1] V. I. ScHAEFER, Chem. Rev. 44, 291 (1949). 
[2] B. Vonnecut, Chem. Rev. 44, 277 (1949). 
H. R. Pruppacuer and R. SAncER, ZAMP 6, 407 (1955). 
[4] H. K. Wercxmann, Arch. Meteor. Geophys. Biokl. [A] 4, 309 (1951). 
[5] L. Rover, Bull. Soc. franc. Minéral. 57, 7 (1928). 
[6] G. FrrEDEL, Legons de cristallographie (Berger Levraut, 1926), 519 ff. 
[7] W. A. WEYL, J. Colloid Sci. 6, 389 (1951). 
[8] S. S. Brrstern and C. E. ANDERSON, J. Meteor. 72, 68 (1955). 
[9] C. L. Hosier, J. Meteor, 8, 326 (1951). 
[10] R. Montmory, Bull. Obs. Puy de Déme 4, 108 (1955). 
Strukturbericht, Bd. II (Z. Kristallogr. 1928-1932) (Akademische Verlags- 
gesellschaft, Leipzig 1937), p. 11. 
[12] H. R. PruppacHer and R. SANGER, ZAMP 6, 485 (1955). 


Zusammenfassung 


Die Experimente fiihren zur Annahme, dass die Eiskeimfahigkeit von den 
Oberflachenkraften herriihrt. Die Epitaxie konnte das dusserst rasche Vereisen 
unterkiihlter Troépfchen nur in extrem giinstigen Fallen erklaren, die aber bei 
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verschiedenen stark keimfahigen Substanzen gar nicht zutreffen. Deshalb wird 
versucht, zusatzlich die Weylsche Theorie der elektrischen Oberflachendoppel- 
schicht heranzuziehen. Tatsachlich weisen die aktiven Stoffe den zur Doppel- 
schichtbildung nétigen starken Unterschied in der Polarisierbarkeit von Anion 
und Kation auf. Die Doppelschicht der Keimsubstanz zerstort diejenige der 
Wasseroberflache, welche die Gitterbildung (Vereisung) verhinderte. 


(Received: March 12, 1956.) 


Uber Verfahren von Steffensen und Householder zur 
Konvergenzverbesserung von Iterationen 


Mauro PICONE zum 70. Geburtstag 


Von ALEXANDER OsTROWSKI, Basel?) 


1. Um die Fixpunkte einer reellen Transformation y = g(x), das heisst die 
Gréssen ¢ mit € = g(¢) zu bestimmen, wird oft das sogenannte Iterationsver- 
fahren?) benutzt, das, mit einem Anfangswert x) beginnend, vermége der 
Rekursionsvorschrift x, ,, = g(*,) in giinstigen Fallen eine unendliche Folge x, 
liefert, die gegen ¢ konvergiert. Gibt es eine Umgebung von ¢ derart, dass fiir — 
jedes x, aus dieser Umgebung die zugehorige Folge x, existiert und gegen € 
konvergiert, so heisst nach Ritt € ein Anziehungspunkt fiir die obige Iteration. 
Existiert dagegen eine Umgebung von ¢ derart, dass fiir jedes x) + ¢ aus dieser 
Umgebung die zugehdrige unendliche Folge x, entweder nicht existiert oder, 
sofern keines der x, den Wert ¢ annimmt, nicht gegen ¢ konvergiert, so heisst 
nach Ritt ¢ ein Abstossungspunkt fiir die obige Iteration. 


2. Existiert ’(¢), so ist ¢ ein Anziehungspunkt, wenn | g’(¢) | < 1 ist, und 
ein Abstossungspunkt, wenn | q’(¢)| > 1 ist. Insbesondere ist im Falle 
| (2) | < 1 die Konvergenz um so schneller, je kleiner | y’(¢) | ist, wenn auch, 
solange g’(¢) + 0 ist, nur «linear». Ist aber g’(¢) = 0, so ist die Konvergenz — 
wesentlich schneller, «iiberlinear», und wenn g(x) — €= O(|x—¢ |”), p>1, ist, 
sogar wenigstens von der Ordnung p, also zum Beispiel fiir = 2 «quadratisch». 


3. Um die Konvergenzverhaltnisse in der Umgebung von ¢ zu verbessern, 
hat STEFFENSEN 1933%) folgendes Verfahren vorgeschlagen: Man bilde von 
Xq = Yo ausgehend die Werte (x) = x, y(x,) = x. und wende sodann die 


1) Mathematische Anstalt der Universitat Basel. 


2) Vergleiche zum Beispiel Fr. A. Witers, Methoden der praktischen Analysis, 2. Aufl. (Wal- 
ter de Gruyter, Berlin 1950), S. 256. 


3) J. F. STEFFENSEN, Skand. Aktuarietidskr. 16, 64-72 (1938). 
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Regula falsi auf die Gleichung F(x) = x — g(x) ~ 0 und auf die beiden Appro- 
ximationen %», x, an. Dann erhalten wir als die nachste Approximation 


%q F(a) — #, F(4%) 


Xy ¥y — 47 


at F(x) — (25) i Hy.— 2 Hy nF , 


Driicken wir hier y, direkt durch yp = x aus, so erhalten wir Vy, = Dy), wo 


VP(P(Y)) = 
¥ — 29(¥) + (Hy) (1) 


ist, wobei allerdings ®(y) = ¢ wird, wenn g(x) linear, = ax + b ist. 

STEFFENSEN betrachtete verschiedene Beispiele, in denen die Iteration ver- 
moége der Funktion @ sehr schnell gegen ¢ konvergiert, sogar dann, wenn 
| p’(f) | = List. WiLLERs bemerkt in seinem Lehrbuch‘), dass «dieses Verfahren 
immer zum Ziel fiihrt», allerdings ohne diese Bemerkung zu begriinden oder 
ihren Giiltigkeitsbereich zu charakterisieren, und gibt in einem Aufsatz®) 
eine elegante geometrische Diskussion des Verfahrens. 


4. Im folgenden werden wir nun beweisen, dass unter sehr allgemeinen 
Regularitatsbedingungen fiir g(x), die wohl praktisch zumeist erfiillt sind, in 
der Tat jeder Fixpunkt ¢ von g(x) ein Anziehungspunkt fiir die Iteration 
vermége der Iterationsfunktion ®(y) ist. 

Es wird dabei im folgenden durchweg vorausgesetzt, dass alle vorkommen- 
den Funktionen und Variablen reell sind. 

Es ergibt sich aus unsern Resultaten, dass, wenn g’(C) existiert und + 1 ist, 
€ ein Anziehungspunkt fiir O(y) mit iiberlinearer Konvergenz ist. Ist anderer- 
seits y’(€) = 1 und gilt zugleich fiir ein A > 1 fir x | ¢ und x} € (x >€ + 0) 


lols) —1| weg lati (et + 0), (2) 


so ist £ wiederum ein Anziehungspunkt fiir die obige Iteration, und es gilt 
. 1 
Die to (3) 


so dass die Konvergenz der neuen Iteration linear ist. Dariiber hinaus ergibt 
sich aber unter geeigneten Annahmen eine weitere Verbesserung des Konver- 
genztypus auch dann, wenn die Konvergenz fir die Iteration vermége v(x) von 
_ der Ordnung f > 1 ist. 


4) Fr. A. WittERS, Methoden der praktischen Analysis, 2. Aufl. (Walter de Gruyter, Berlin 


1950), S. 259. 
5) Fr. A. Witters, Z. angew. Math. Mech. 28, 125—126 (1948). 
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5. In seinem 1953 erschienenen Buche®) hat A. S. HOUSEHOLDER das Stef- 
fensensche Verfahren dahin verallgemeinert, dass er zugleich zwer Iterationen 


Y= Q(x), Y= Pe(X) 


mit gemeinsamem Fixpunkt ¢ betrachtet und aus ihnen die neue Iterations- 


funktion 
Bly) = —2 Pill) — Paty 


(VY) Po(¥) 
) PalY) 4 
y— P(¥) — aly) +9 a (y)) ( ) 


bildet. Doch bezieht sich seine Diskussion nur auf den Fall, wo 


(pi(2) — 1) (p22) — 1) + 0 (5) 


ist, und y,(x) und ,(x) geniigend oft differenzierbar sind. 

Im folgenden behandeln wir nun gleich den allgemeineren Householderschen 
Ansatz unter wesentlich weniger einschrankenden Voraussetzungen. Es ergibt 
sich dabei allerdings, dass die Benutzung von zwei verschiedenen Iterations- 


funktionen @,, y, nur in einem recht speziellen Fall (Nr. 12) ein besseres Er- — 


gebnis liefert, als wenn man eine geeignete dieser Funktionen direkt im Ansatz 
von STEFFENSEN benutzt. Die genauen Formulierungen unserer Resultate 
findet man weiter unten in den Nummern 7, 12-16, 20. 


Es sei endlich noch bemerkt, dass wenn @j(C) und (C) beide positiv sind, 


man sich in den weiteren Betrachtungen sowohl in den Voraussetzungen als 
auch in den Behauptungen entweder durchweg auf die rechtsseitige oder durch- 
weg auf die linksseitige Konvergenz von y gegen ¢ beschranken kann. 


6. Wir bemerken zuerst, dass die Bildung von ®(y) vermége (4) invariant 
ist gegentiber der Translation, so dass man bei den weiteren Beweisen ¢ = 0 
annehmen darf. In der Tat verifiziert man unmittelbar, dass, wenn man 


pilx) = (C+ x)—¢ (= 1,2) 


setzt und die zu den beiden Funktionen y,(x) und ,(«) gehérende Househol- 
dersche Iterationsfunktion mit Y bezeichnet, dann die Relation gilt 
(ys Dy Pili 
Aus (4) folgt sofort: 
Oly) — y = —— Al = SD 
"= ul) — al) + PilGalV)) o 


e ep en HousEHoLpeEr, Principles of Numerical Analysis (McGraw-Hill, New York 1953), 
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i. Da im folgenden durchweg die Existenz der ersten Ableitungen von 
yi(y) in ¢ vorausgesetzt wird, gelten also fiir y > C, |v —¢| = @ gesetzt, die 
Formeln 


Pil¥) = 4 (y—C) +0(e) = 1,2; y>2), (7) 


wo «; die Ableitungen von @,(y) in € sind. 
Gilt nun (5), das heisst 


(a, — 1) (%—1) +0, (8) 


so lisst sich bereits beweisen, dass fiir die Iteration vermige P(y), C einen An- 
ziehungspunkt mit tiberlinearer Konvergenz darstellt. In der Tat ergibt sich aus 
(4) und (7) fiir ¢= 0: 


P(y) Oy Oty YR + 0(Q*) — (8, 9 0(@)) (> eae 0(@)) bz 
Y— & V— A V+ A & Y + O(0) 
0(0°) 


(I — a) (i s«,) y+ of) o(@) - 


8. Will man allerdings einerseits die Falle «, = 1, «, = 1 betrachten, und 
andererseits auch in dem oben zuerst betrachteten Fall weitere Aussagen 
machen, so miissen tiber das Verhalten von @,(y) weitere Annahmen gemacht 
werden. Wir formulieren diese Annahmen folgendermassen: 

Es soll fiir 1 = 1, 2 je eine Zahl A; > 1 existieren derart, dass fiir y >¢ 


Ey) = #40 — 9-9) oa) a, >1Ly>ti=1,2) (9) 
a) 7 


gilt, das heisst 
gly) = a(y — 6) + Ely) 0”, (10) 


wo E,(y) fiir y > C beschréinkt bleibt. 


Wir nehmen in der folgenden Diskussion an, dass A, und A, endlich sind, das 
heisst, dass die y,(y) nicht linear sind. Auf den Fall, dass p, oder @, linear sind, 
ist es offenbar nicht notig einzugehen, da dann ¢ aus einer linearen Gleichung 
zu bestimmen ist. 5 


9. Wir bezeichnen nun den Zahler und Nenner im Ausdruck (4) bzw. mit 
Z(y) und N(y), so dass 


Z(y) = ¥ G1 (Pal¥)) — Pry) Pal) » (11) 


N(y) = ¥ — galv) — Holy) + Pr (G2) (12) 


ist. Ferner soll im folgenden von der Konvention Gebrauch gemacht werden, 


AY 
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dass, wenn bei den Funktionsbezeichnungen E,, 9; kein Argument angegeben 
wird, dieses Argument immer y ist. 


10. Setzen wir in Z(y) die Ausdriicke (10) ein, so ergibt sich fiir @ = 0 
Zy) = ety tty Y? + oy Epo’ y + Ey (ps) | Go|" ¥ — 4 He ¥? 
= a y Eye” — afy Eye" = £, Eg = 
= E, (92) | P2|" ¥ — % y By 9" — E, Ey". 
Hier gilt aber | y, | = | % |e + O(0”), 
- | Pal" = [arg|* g + O(Q"**) , (13) 
so dass wir schliesslich erhalten 
Z(y) = T, yo" + O(o"**) , (14) 
T, = |a&9{" Ey (Qo) — a Ey. (15) 


Setzen wir andererseits (10) in den obigen Ausdruck (12) fiir N(¥) ein, so 
ergibt sich fiir @ = 0 ; 


N(y) = (1 — 4) (1 — a) y — Ey 0" — Ey o* + a Ey 0 + Ey (2) |p|". (16) 9 
Machen wir hier von (13) Gebrauch, so erhalten wir 
N(y) = (64 — 1) (ae = 1) y+ (oq —"1) Bo 4D, 0 O02 a 
Ty = |o2|"* E1(p2) — Ey. (18) 
11. Diese Formeln miissen indessen in dem Fall «, = 1, a, = -+ 1 noch 
weiter entwickelt werden. Man beachte, dass fiir «= -- 1 aus (7) folgt, wenn 
sgn y = e gesetzt wird, 
[Pol =o (+ ea, 0%") (¢,=-+1). 
Setzt man dies in (16) ein und beriicksichtigt «, = 1, so ergibt sich 
N = 0" [Ey (gq) (1 + € a Ey oe) — Ey), 
SANT oe 
@ “N= E,(p,) — Ey + €A, a, Ey (2) Ey ot Oe i cee |%2|= 1). (19) 


Da wir bei der Betrachtung der Falle mit «, = | x, | = 1 die weitere Annahme 
machen werden, dass die Ableitung von E,(y) in den beiden einseitigen offenen 
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Umgebungen von ¢ existiert, kbnnen wir hier weiter entwickeln: 
E,(P2) = Ey(a, y) + Eye Ex(é), E=o y+ OE, 0", 0<0<1 ' 


E,( a) — Ey (a4 2 Ey the , 
er aE et EEYE) = O(6 EY)), 


da fiir 9 > 0 zugleich |&| ~ @ ist. Daher folgt wegen (19) 
Ey(P2) = Ey(a, y) + op * O(E Ey(8) , (19a) 
oN = Ey(ay y) — Ey +A, a E,(ay y) Es 0” | 
+ OE EO) +00) (=le|=1), | 
wo € > 0 mit 0 + 0 gilt. 


12. Machen wir nun die Annahme (8), so folgt wegen (9) und (10) aus (17), 
wenn Min (A,, A,) =A > 1 gesetzt wird, 


Nunmehr folgt aber wegen (14) 


Z iit 
as Be te hfe =F ot One) : (21) 

Ist a, = 0, Tz = 0, so folgt aus (14) bereits Z/N = O(o9"**~1). Wir sehen 
daher: 

Unter den Annahmen (8), (9) und (10) gilt B(y) — € = O(0") und fiir a, =0 
sogar D(y) — € = O(o%**-*) . 

Man sieht insbesondere, dass fiir «, + 0 im Falle (8) es auf die Wahl von 
(2 gar nicht ankommt, so dass man ebensogut y, = @, setzen kann. Ist dagegen 
% = 0, so kann es vorteilhaft sein, y, verschieden von g, zu wahlen, wenn 
namlich A, > A, ist’). 


13. Wir betrachten nunmehr die Falle, in denen (8) nicht erfiillt ist. In 
diesen Fallen miissen zusatzliche Annahmen iiber die Gréssen E,(y) gemacht 
werden. Wir betrachten zunachst den Fall «, + 1, a, = 1 und beweisen: 


7) Die Annahmen unseres oben hergeleiteten Resultats treffen zu auf den von HousEHOLDER 
(siehe Fussnote 6) allein betrachteten Fall. Nimmt man der einfacheren Schreibweise wegen ¢ = 0 
an, so setzt HousEHoLpER die Funktionen ,(*) (¢ = 1, 2) in der Form 


CAA ae de awe dalle ste 


mit gemeinsamem natiirlichem 7 an, wobei fiiry = 1 vorausgesetzt wird, dass (ox?) —1) (0?) —1)+0 
ist. Er erhalt dann fiir > 1: D(y) = 0(07"~') und fiiry =1: W(y) = O(0?). 
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Es sei 0, + 1,4 = 1; es set ferner die Ableitung E,(y) sowohl 1m der rechts- 


seitigen als auch in der linksseitigen Umgebung von C, mit eventueller Ausnahme — 


von €, vorhanden und dort 0 E,(y) gleichmdssig beschrainkt. Ist zugleich mit 
E,(y) auch 1/E,(y) in einer Umgebung von C gleichmdssig beschrankt, so gilt 


Py) — £ = O(0"). ; 
In der Tat gilt nach den obigen Annahmen fiir ¢ = 0 nach dem Mittelwert- 


satz der Differentialrechnung, wegen (15) und (18): 
T, = Ty = Ely + Eso") — Ei = Ole). 
Daher folgt aus (17) 
N = (a, — 1) Eo 4 O(e**4~*) 
und daher schliesslich, wegen (14), 


Z 1 O(9% +42) 3 | 
ees wae = O(g* 22 
2 N (a, — 1) BE, o” 4 O(o%*4-1) (e )4 ‘ ) 


was zu bewelsen war. 
Wir sehen, dass sich fiir a, + 0 eine wesentliche Konvergenzverbesserung 
ergibt, wahrend im Falle «, = 0 wir die gleiche Gréssenordnung fiir ® wie fiir 


y, erhalten haben, und diese Gréssenordnung lasst sich, wie eine eingehendere. : 


Diskussion zeigt, im allgemeinen nicht verbessern. 


14. Wir betrachten nunmehr die Falle, in denen «, = 1 ist. Hier werden 
wir je nach dem Wert von «, mehrere Unterfalle unterscheiden. 
Wir beweisen zuerst: 


Ist a =1, a =0 und 1/E,(y) in einer Umgebung von € gleichmdssig be- 


schrankt, so gilt B(y) — ¢ = 0(@*). 
In der Tat folgt unter den jetzigen Annahmen fiir ¢ = 0 aus (14), (15), (17) 
und (18) 
Z=O(gh**), NS = Bo OG se 


woraus die Behauptung unmittelbar folgt. 


Wie man sieht, ergibt sich auch in diesem Falle keine Verbesserung gegen- — 


liber po. 


15. Ist % = 1, % + 0, + 1 und strebt E,(y) fiir y | € bzw. yt € gegen von 0 
verschiedene Grenzwerte P,, P_, so gilt 


By) —C=Py tole), pa tgalta ae (23) 


eager 


ee 
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In der Tat ergibt sich jetzt fiir € = 0 aus (14), (15), (17) und (18) 


Z~ (cae E4(P) — ay E,) oe o” > P Ate, |” — %) y 0” , 
N~ (|a.|" Ey (2) — Ey) 0 a (|_| —o) o” , 


je nachdem ob y | ¢ oder y } € gilt, und daher Z/N ~ B y, wie behauptet. 

Offenbar ist in unserem Falle 6 + 0, so dass die Kombination von , mit 
@2 auf jeden Fall ein wesentlich schlechteres Ergebnis liefert, als wenn man 
Pz mit gp, kombiniert. Ist x. << —1, so ist B > 1, so dass ¢ fiir ® zu einem Ab- 
stossungspunkt wird. Liegt «, zwischen —1 und 0, so ist 6 negativ und kann 
sowohl < —1 als auch > —1 sein. Immerhin ergibt sich fiir hinreichend 
kleine Werte von a: |6|<1. Ist andererseits « >1, so ist B positiv und kleiner 
als 1, so dass sich immerhin ein Anziehungspunkt ergibt. Ist aber 0<a,<1, 
so gilt sogar 0 < fy < a, so dass eine Konvergenzverbesserung gegeniiber y, 
zu konstatieren ist. In der Tat gilt ja in diesem Falle 


B Od et fl yl A Pa 


My ie | oc, |” 


16. Wir betrachten nunmehr den Fall «, = «,=1. Hier ergibt sich fiir 
€ = 0 aus (20), wenn £j(y) in den beiden einseitigen offenen Umgebungen von 
€ existiert und 9 F{(y) mit @ > 0 gegen 0 konvergiert, 


o-* N = eh E, E, 0% [1 + o(1)]. 
Daher folgt aus (6) und (10) 


EE outs y 
Do = oot = —_ + 0(0) . 
f CAPE pepo rT hs ol} Ay (2) 


Wir sehen: 
Tst = %=1, E,E,+0 und ist @ Ej(y) +0 fiir € >0, so gilt: 


-C=(1-5) -tey-o. (24) 


In diesem Falle erhalten wir also die lineare Konvergenz. Uber die Bedeutung 
der Annahme 0 Ej(y) > 0 vergleiche man Nr. 23. 


17. Endlich sei der Fall «, = 1, «, = —1 betrachtet. Hier folgt aus (20), wenn 
wir annehmen, dass E/ vorhanden ist und 9 E;(y) = O(1) bleibt fiir 0 > 0: 


o-* N = E, (—y) — E,(y) + Oe") - (25) 


ZAMP VII/15 
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Andererseits folgt aus (14) und (15) durch Division mit y o”: 


mie ' 
re = E,(p2) + Ex(y) + O(o” ’) 


und daher, wegen der Beschranktheit von @ Ej, nach (19a) 


mi, 
Big) 


yoo E,(— y) + Ey(y) + Ofe"*>}).. (26) 


Fiir die weitere Diskussion ist offenbar das Verhalten von £,(—¥y) und 
E,(y) massgebend. Hier betrachten wir zwei Falle, die wir als den Fall u) und 
Fall g) bezeichnen, weil sie im Falle eines ganzzahligen A, mit der Annahme 
zusammenhangen, dass A, ungerade bzw. gerade ist. 


18. Im Falle ) setzen wir voraus, dass fiir y > 0 die Relation besteht 
u) E,(y)=eP+O(o) (P +09), ; (27) 


€ = sgny gesetzt. In diesem Falle folgt aus (25) und (26) 


“= = 0(0) + Oe"), bea" NS ee Pa Ole) Ces: 
und daher durch Division 


u) B= X = O(0%) + O(o%) . (28) 


19. Im Falle g) dagegen setzen wir voraus, dass 


8)» Exly= P Ole) “(Pe 0), (29) 
ist. Hier folgt 


i = Sy 
eg * N= Ole) + Oe"), f= = 2P + O(e) + O(e*"), 


und daher durch Division 


®D  2P+O0(0) + Ofo—}) 
y O(a) + O(e*-1) 


(30) 
In diesem Falle gilt | B/y | +00, so dass ¢ sicherlich kein Anziehungspunkt 
fiir P(y) ist. Man kann sogar im Fall g) Beispiele bilden, in denen ¢ nicht einmal — 
ein Fixpunkt bleibt. Setzt man zum Beispiel Q(x) = 4+ x?, (x) =— 4% + 28, . 
so ergibt eine leichte Rechnung, dass ® = y — 1/y ist, so dass ® fir y >0 | 
nicht endlich bleibt, anstatt gegen 0 zu konvergieren. 


SS 
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20. Zusammenfassend erhalten wir: 

Es set = —a@=1 und es sei Ej(y) vorhanden und 0 Ex(y) gleichmdssig 
beschrankt, sowohl fiir y+ € als auch fiir y | ¢. Nimmt man an, dass fiir ein 
P+0 E,(y) = P+ O(@) oder = —P + O(o) ist, je nachdem, ob y von rechts 
oder von links gegen ¢ konvergiert, so gilt (28), so dass die Ordnung der Konvergenz 
der Iteration vermige P wenigstens gleich Min (2, Ay) > 1 ist. 

Gilt dagegen fiir beidseitige Konvergenz von y gegen C, E,(y) = P+ O(o) fiir 
ein P+ 0, so ist € sicherlich kein Anziehungspunkt fiir die Iteration vermige D 
und braucht nicht einmal ein Fixpunkt dieser Iteration zu sein. 


21. Man wird sich nun vor allem dafiir interessieren, ob man durch den 
Householderschen Ansatz ein besseres Ergebnis erhalten kann, als indem man 
eine geeignete der beiden Funktionen 9, g, mit sich selbst, wie im Ansatz von 
STEFFENSEN, kombiniert. Geht man die oben betrachteten Falle durch, so sieht 
man, dass in allen Fallen, bis auf den in Nr. 12 betrachteten Fall, man durch 
den Steffensenschen Ansatz entweder ein besseres Resultat bekommt oder ein 
ebensogutes mit weniger Annahmen. Eine Verbesserung gegentiber dem Stef- 
fensenschen Ansatz ist im Falle der Nr, 12 dann vorhanden, wenn A, >A, und 
a = 0, a, + O ist, da dann beim Householderschen Ansatz die Ordnung der 
Konvergenz A, + A, — 1 ist, wahrend sie, wenn g, mit sich selbst kombiniert 
wird, nur A,, und wenn q@, mit sich selbst kombiniert wird, nur 2 A, — 1 ist. 


22. Es mége noch der Sinn der in der obigen Diskussion benutzten Annah- 
men iiber Ej(y) erlautert werden. Die Annahme, dass @ Fj(y) = O(1) ist, lauft 
darauf hinaus, dass i(y) in einer Umgebung von C vorhanden und dass dort 


p(y) — % =O(o""") (y+) (31) 


ist. In der Tat folgt aus (10), dass E,(y) fiir y + ¢ zugleich mit ¢,(y) in einer 
Umgebung von ¢ differenzierbar ist. Differenziert man nun (10) nach y, so 
folgt fiir ¢ = 0 


gi(y) =o + Ey(y) 0% + 4, Ey(y) oe * (e= sgn y,, (32) 


und daher nach (31) 
e E,(y) = O(1) . 


Existiert umgekehrt E{(y) und ist @ E{(y) = O(1), so folgt aus (32) die Re- 
lation (31). 


23. Was andererseits die Bedingung 
0 Ex(y) >0 (9 +9) (33) 
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anbetrifft, so ist sie etwas allgemeiner als die Annahme, dass 


noe? aces 
ae 
|y_+0 rt ope 


Man kann namlich leicht zeigen, dass aus (33) und der oben vorausgesetzten 


Beschranktheit von E,(y) nicht etwa folgt, dass E,(y) fir y }¢, yt ¢ gegen 
von 0 verschiedene Grenzwerte konvergiert. Wird aber neben (33) auch noch 


explizite vorausgesetzt : 


(34) 


| o, #0 (y{¢), 


o> (35) 
|-c_ +0 (yt ¢), 


E,(y) 


so lasst sich zeigen, dass die Annahmen (33) und (35) zusammen 4quivalent 
sind mit der Annahme (34), wobei 


74 =A 4. _ (36) 


In der Tat folgt aus (32), wenn man dort (33) und (35) beriicksichtigt, sofort 
(34) und (36). 
Gilt aber umgekehrt (34), so wende man auf den Quotienten 


E . 
tee 


die Bernoulli-l’Hépitalsche Regel an; dann folgt aus (34), wegen (E,(y) ot)’ = 
= ily) — % 
E,(y) — oe ’ 
1 
das heisst (35) und (36). Wird aber dies in (32) eingetragen, so ergibt sich 
sodann auch (33). 


24. Das Steffensensche Verfahren wird gelegentlich in der neueren Literatur 
als «Aitken’s 6?-Prozess» bezeichnet, was allerdings, wie uns scheint, auf einem 
Missverstandnis beruht. In der betreffenden Arbeit von AITKEN®) handelt es 
sich um die Bemerkung, dass, wenn die Zahlenfolge a, aus einer linearen Rekur- 
sionsvorschrift endlicher Ordnung (einer linearen Differenzengleichung mit 
konstanten Koeffizienten) entspringt, und a, 41/4, >A gilt, dann der Ausdruck 


2 
G41 U1 — UW 
Ay = 2 ay + ay_1 


(37) 


8 aes F 
: ) A. C. AITKEN, Studies im Practical Mathematics. 11. The Evaluation of the Latent Roots and 
atent Vectors of a Matrix, Proc. roy. Soc., Edinburgh, 57, 291-295 (1937). 


we 
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noch schneller gegen A konvergiert. Dies trifft in der Tat unter sehr allgemeinen 
Annahmen zu, wobei indessen festzustellen ist, dass, wenn die Konvergenz von 
a, .1/a,, wie es gewohnlich der Fall ist, linear ist, auch der Ausdruck (37) gleich- 
falls linear gegen A konvergiert, aber mit kleinerem Quotienten als a,,,/a,. Es 
wird also dabei die Ordnung der Konvergenz nicht erhoht, und andererseits wird 
auch durch die Anwendung des Aitkenschen 6?-Prozesses die ganze Folge a, 
nicht gedndert, sondern nur in vorteilhafterer Weise ausgeniitzt. Wird anderer- 
seits durch den Steffensenschen Iterationsprozess die Iterationsfunktion (x) 
durch die Iterationsfunktion ®(y) ersetzt, so wird dabei auch die ganze Folge 
der Iterationswerte vollstandig geaindert. Insofern scheint es uns, dass, ob- 
gleich die Formel (37) formale Ahnlichkeit mit der Formel (1) aufweist, es 
kaum gerechtfertigt sein diirfte, den Steffensenschen Ansatz mit dem Aitken- 
schen 6-Prozess zu identifizieren. 


Summary 


STEFFENSEN discovered that the iteration by y = (¥%) can be considerably 
improved if g(*) is replaced by the function 


y p(p(y)) — ply)? 
PY) = Fa eint oP) 


This procedure has been generalized by HousSEHOLDER, who starts from two 
iterations y = 9,(*), y = 9,(*) with a common fix-point and builds up the new 
iterating function 


Lean! 91(H2(¥)) — 91(V) P2(¥V) 
gre PY) — ald) + r(H2()) © 


We prove under very general conditions that any fix-point of (*), resp. any 
common fix-point of g,(¥) and g,(*) becomes a point of attraction for Stef- 
fensen’s or Householder’s iterating function, and compare both procedures. 


(Eingegangen: 26. August 1955.) 
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Thermal Stresses in Viscoelastic Structures’) 
By Witt1aM Pracer, Providence, Rhode Island, USA.?) 


1. Introduction 


In the standard texts on the theory of structures, more space is allotted to 
trusses than seems warranted by their role in modern structural engineering. 
This is done because the discussion of the mechanical behavior of trusses 
requires but a minimum of mathematics and thus provides an excellent op- 
portunity for developing the student’s understanding of the mechanical be~ 
havior of structures in general. In a similar manner, the influence of various 
mechanical effects on structural behavior is most readily. investigated for 
trusses. As a rule, the extension of the results of such an investigation to other 
types of structures does not present any difficulties. 

In recent papers FREUDENTHAL [1], [2]%) used the simple truss shown in 
Figure 1 as a convenient structural model for the discussion of inelastic thermal 


A B C 


/ 


(e) 


Pp 
Figure 1 
FREUDENTHAL’S structural model. 


stresses. The state of loading of this simply indeterminate structure was sup- 
posed to be specified by a single parameter, the load intensity P, and only the 
bar OB was supposed to exhibit inelastic behavior; the two other bars were 
assumed to be perfectly elastic. 

Even so simple a model may reflect many characteristic traits, but it cannot 
be expected to reveal all these traits*). A more general investigation of thermal. 
stresses in inelastic trusses seems therefore desirable. 


1 : . . . 
) The results presented in this paper were obtained in the course of research sponsored by the 


Ballistic Research Laboratories of Aberdeen Proving Ground under Contract DA-19-020-ORD-798. 
2) Brown University. 


8) Numbers in brackets refer to References, page 237. 


4 ; P 
) For perfectly plastic structures, for instance, this intuitive statement is confirmed by the 
work of PRAGER and Symonps [3], [4]. 


Vol. VII, 1956 Thermal Stresses in Viscoelastic Structures 231 


The present paper is concerned with the thermo-mechanical behavior of a 
statically indeterminate truss that consists of Maxwell bars. The effects of 
typical variations of loads and temperatures on the stresses in the truss are 
discussed. While the actual relation between creep rate and stress in structural 
metals is less simple than the linear relation stipulated for a Maxwell material, 
the structural model considered in the following may be expected to indicate 
the general effect of creep on the stresses in indeterminate structures. 


2. Fundamental Relations 


The principle of virtual displacements is the only tool of structural theory 
used in the following discussion. For a truss, this principle involves, on the 
static side, loads P and bar forces S that are in equilibrium with these loads, 
and, on the kinematic side, joint displacements u and bar elongations A that 
result from these displacements. The principle of virtual displacements then 


asserts that 
DJA Pet, (1) 


where the sum on the left includes all bars and that on the right all joints of 
the truss, the dot being used to indicate the scalar product. 

The discussion of the mechanical behavior of highly indeterminate structures 
is greatly facilitated by the introduction of orthonormal states of stress (see, for 
instance, [5]). Consider two states of stress in a given truss, and let S’ and S” 
denote the corresponding forces in a generic bar. If / denotes the length of this 
bar, A its cross-sectional area, E its modulus of elasticity, and c = //(AE) its 
elastic compliance, the two states of stress are said to be orthogonal if 


Sra SeGl= 0% (2) 


where the sum includes all bars of the truss. If moreover 


DP? GS Paid ON ge (3) 


the two states of stress are said to be orthonormal. Two groups of loads acting 
on a truss are called orthonormal, if the elastic bar forces corresponding to them 
are orthonormal. These elastic forces must be computed by considering only 
the elasticity of the bars but neglecting all other rheological properties, e.g. 
viscosity or plasticity. 

A system of loads that depends on m parameters is most conveniently 
treated as a linear combination of m orthonormal groups of loads. The coef- 
ficient p/ of the i-th group will be called the intensity of this group in the 
linear combination. If S/ is the elastic force caused in a generic bar by the 
unit intensity of the i-th group, the elastic response of this bar to an arbitrary 
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state of loading of the considered type is given by 
SO =p Si + Pe Sat ++ + bm Sm: (4) 


This elastic response, of course, represents the actual force in the bar only if 
all bars of the truss are perfectly elastic and thermal stresses are absent. 

When inelastic behavior or thermal stresses are to be considered, the 
deviation of the actual state of stress from the elastic response to the given 
loads constitutes a state of residual stress®). In a truss with m redundancies, 
the most general state of residual stress can be treated as a linear combination ~ 
of m orthonormal states of residual stress. Let p;’ be the intensity of the j-th of 
these states in the linear combination and Sj’ the residual force caused in a 
generic bar by the unit intensity of the j-th state of residual stress. The actual 
force in this bar can thus be written in the form 


Sap Sites 4 pl Si + pl SE +--+ + PF SK. (5) 


It is readily shown that any elastic state of stress is orthogonal to any 
residual state of stress. Indeed, the joint displacements w’ associated with an 
elastic state of stress S’ produce the bar elongations 4’ = c S’, and the bar — 
forces S” of a residual state of stress S” do not require any loads P” for 
equilibrium. Application of the principle of virtual displacements to the kine- 
matic variables u’ and 4’ =c S’, and the static variables P” =0 and S” — 
therefore shows the states S’ and S” to be orthogonal. It follows from this ~ 
that the m+ n states involved in the linear combination (5) form an ortho- 
normal set. 


3. Thermal Stresses in Elastic Trusses 


At the time ¢ = 0, let the truss be free from loads and stresses. The absence 
of loads implies that the intensities #1, 3, ... , £), of the m orthonormal load 
groups vanish initially, but the absence of thermal stresses does not imply 
that all bars have the same initial temperature. Any subsequent thermal 
stresses, however, will depend only on the differences 6 between the tempera- — 
tures of the various bars in the considered state and the initial state. For 
t > 0, these temperature differences and the intensities of the m orthonormal 
load groups are supposed to be given as functions of time; these functions are 
supposed to vary slowly with time, so that a quasi-static analysis is appropriate. 

Let S be the force in a generic bar and 6 the difference between the 
temperatures of this bar in the instantaneous and initial states. The total 
elongation A of the bar now consists of the elastic part 4 = c S and the thermal 
part 2 = «16, where « is the coefficient of thermal expansion. Since the 


5 : : ; 
) This term is here used as the equivalent of the German “Selbstspannungszustand’. 
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total elongation A= c S + «16 of the bars correspond to a system of joint 
displacements, the principle of virtual displacements now furnishes 


D> (¢S+al6)S"=0 G =T,2,...,n). (6) 


Here S is given by (5). The states appearing in (5) are orthonormal; substitution 
of (5) into (6) therefore furnishes 


bj =—D7 al S! C= Fr 2m BM) (7) 
In the expression (5) for a generic bar force, the intensities #/ (i = 1, 2, ... , m) 
are given as functions of time, whereas the intensities Dit se 2a, a) 


must be found from the given temperature differences 6(t) by means of (7). 

It is convenient to introduce the fictitious thermal forces @ that the 
temperature changes @ would produce if the bars were not allowed to extend 
or contract. Since 

ga (8) 


c 
equation (7) can be written in the form 
pp =) OS; sean Wi AO (9) 


where the superscript (e) has been used to indicate that (9) refers to a perfectly 
elastic truss. 

Equation (9) shows that there will not be any thermal stresses, if the 
fictitious state of stress specified by the forces (8) is orthogonal to all residual 
states of stress S/’. 


4. Thermal Stresses in Homogeneous Trusses with Maxwell Bars 


The rate of elongation of a Maxwell bar is given by 
. . ‘Si . 
A=e(S+2-6), (10) 


where the dot indicates differentiation with respect to time and S is the force 
in the bar, 7 the relaxation time treated as independent of temperature and 
stress, and @ the fictitious force (8). 

Since the rates of elongation of the bars of the truss correspond to a system 
of joint velocities, the principle of virtual velocities furnishes 


Ye(S+2-6)sf=0 @=1,2,...,9). (11) 


In discussing equation (11), assume first that the relaxation time t has the 
same value for each bar (homogeneous truss). Substitution of (5) into (11) and 
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use of the orthonormality of the states appearing in (5) then leads to following 


differential equation: 
pie = = pf (te; Doe NS (12) 


where p/’ is given by (9). The solution of (12) that satisfies the initial con- 
dition ~/’ = 0 for t = 0 is given by 


t 
pf = {BLOW el at (13) 
0 


When all bars have the same relaxation time (homogeneous truss) there is 


therefore no interaction between the thermal stresses and the stresses caused 
by the external loads. In particular, in the absence of thermal effects the stresses 
that given fixed loads produce in a homogeneous truss with Maxwell bars 
equal the stresses that these loads produce in the corresponding elastic truss. 
The analoguous theorem for the viscoelastic continuum is due to ALFREY [6]. 

Assume that /) varies with time as shown by the line OAB in Figure 2. 
It then follows from equation (13) that #/ reaches its maximum at the instant 
t=" and that 


pi) = a Ae Fh). (14) 


Figure 3 shows #;'(¢*) as a function of the ratio ¢*/r. It is seen that the viscosity 
of the Maxwell bars leads to a substantial attenuation of the thermal stresses 
predicted by elastic theory except when ¢*/r <1, i.e. when the temperature 
changes are effected during an interval of time ¢* that is much smaller than 
the relaxation time T. 


Figure 2 
Assumed variation of elastic thermal stress 


Figure = 
Dependence of maximum actual thermal 
pr) with time ft. stress p/(¢*) on ratio between warm-up time 
t* (see Figure 2) and relaxation time T. 
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Another important type of variation of p; with time is represented by 
the line OAB in Figure 4: p/) is suddenly raised from 0 to 1 at t= 0 and is 
thereafter made to fluctuate according to p/ =1+ asinw t. The correspond- 
ing variation of #/ is found to be 


| RT OP Aw awWTt : 
p; =e (1 cea 4 te ait (cosmt+artsinwt). (15) 


The first term on the right-hand side of (15) represents the transient effects 
of the sudden raise of #/ from 0 to 1 and the start of the fluctuations in 


Figure 4 


Response of the actual thermal stress p} to the fluctuations of elastic thermal stress represented 


by the line OA B. (T = period of elastic thermal stress, T = relaxation time.) 


pi at t= 0; the second term on the right of equation (15) represents the 
steady state that is reached when the transient effects have died out. Whereas 
p! © fluctuates about the value #/ = 1, the steady state response of $/ 
fluctuates about the value #/ = 0. This remark is important when fatigue has 
to be considered. 

Figure 4 shows the variation of #;’ with time for some values of the ratio 
between the relaxation time t and the period T = 2./ of the fluctuations. 
In the limit t/T +00, we have $j’ = p//. The smaller t/T, the faster does 
the steady state become established. In the steady state the ratio between 
the amplitudes of // and p/ is given by wt (1 + w? 1?) 71. 


5. Thermal Stresses in Inhomogeneous Trusses with Maxwell Bars 


Consider next an inhomogeneous truss in which various bars have different 
relaxation times, either because they are made of different materials or because 
their average temperatures differ so much that the dependence of the relaxation 
time on the average temperature must be taken into account. Substitution 
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of (5) into (11) and use of the orthonormality of the states in (5) then yields 


b; ae b,; 2; eee Di b= p ay ay; Py ri: Aiee st Iai Pn | 


(16) 
G22, (oa | 


where aad putz 
elles SSE) ee (eae (17) 


The general solution of the system (16) is the sum of a particular solution 
of this system and the general solution of the corresponding homogeneous 
system. The latter describes the relaxation of a generic state of residual stress 
in the absence of loads and temperature changes. The discussion of this complex 
relaxation process is simplified by the introduction of the characteristic states 
of residual stress. A state of residual stress is called characteristic if, in the 
absence of loads and temperature changes, the forces in all bars decrease as 
e—/*, where the relaxation time t associated with the considered characteristic 
state has the same value for all bars of the truss. It can be shown®) that an 
arbitrary state of residual stress in a truss with m redundancies can be treated 
as a linear combination of n orthonormal characteristic states of residual stress. 
In general, the relaxation times 7}, T., ... , T, associated with the characteristic 
states are distinct, but in special circumstances two or more of these states 
may have the same relaxation times. For example, for the homogeneous truss 
considered in the preceding section, all states of residual stress have the same 
relaxation time T. 

When the » arbitrary orthonormal states of residual stress used in the 
preceding analysis are replaced by orthonormal characteristic states of re- 
sidual stress, equation (16) reduces to 


b+ 2t = f° _ ap —- Sta io GS 1, een 


7 
where 


ps MeO OCIS? ox = Sy (SESE) (19) 


In equations (18) and (19), #;’ denotes the intensity of the i-th characteristic 
state, T; its relaxation time, and S;’ the force in a generic bar caused by the 
unit intensity of the 7-th characteristic state, and the fictitious thermal force 
@ is given by (8). 

The first term on the right-hand side of (18) represents the thermal effects, 


whereas the other terms depend on the external loads. Comparison of (18) 
6 ; : 
) The proof of this statement, contained in the first draft of the present paper, was deleted 


after the author’s attention had been drawn to a i i 
; : : recent paper by Bior [7], w 
discussion of linear relaxation processes. paper by [7], which contains a general 
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and (12) shows that the thermal effects in an inhomogeneous truss resemble 
those in a homogeneous truss except that the various states of residual stress 
may now be associated with different relaxation times. 

To illustrate the behavior predicted by the load terms on the right-hand 
side of equation (18) in the absence Of thermal effects, assume that all loads 
are suddenly brought on the truss and thereafter kept constant. The immediate 
response of the truss to this sudden loading is purely elastic, i.e. immediately 
after the loading #;/=0 (¢=1, 2, ..., ). The intensity #/’ of the i-th 
characteristic state of residual stress then tends towards the asymptotic value 
obtained by multiplying the load terms on the right-hand side of equation (18) 
by t,. If this asymptotic value is assumed to be unity, the variation of #;' 
with time is given by p/ =1—e-‘";. Since the characteristic states with 
small relaxation times approach their asymptotic intensities faster than the 
characteristic states with large relaxation times, the transition from the elastic 
state of stress prevailing immediately after loading to the ultimate state of 
stress is rather complex. A graphical treatment of this transition in a special 
case has been given by MEACHAM [8]. 

An interesting example of combined load and temperature effects is ob- 
tained by assuming the loads to vary with cosw ¢ and the temperature differ- 
ences with sinw ¢. The loads then have maximum intensity when the tempera- 
ture differences vanish and vice versa. The right-hand side of equation (18) is 
of the form A cos ¢ and the steady state behavior of p/’ is given by 


~; = Bcos(wi— 9), (20) 
where 
,A 
B — Mayu nee ) tan = W T; . “ih 
i+ oft P=o (21) 


Note that the amplitude A of the right-hand side of equation (18) reflects 
the temperature fluctuations as well as the load fluctuations. As far as the 
residual stresses are concerned, temperature fluctuations can therefore be re- 
placed by equivalent load fluctuations. Since, however, the intensities bi, Pas 
pi}, -.. , Pi, of the elastic stresses depend only on the loads, such an equivalence 
does not exist for the total stresses. 
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Zusammenfassung 


Das thermo-mechanische Verhalten eines statisch unbestimmten Fachwerks 
wird untersucht unter der Annahme, dass die Fachwerkstabe der Maxwellschen 
Spannungs-Formanderungsbeziehung gehorchen. Insbesondere wird der Einfluss 
typischer Last- und Temperaturanderungen besprochen. Obwohl der genaue Zu- 
sammenhang zwischen Spannung und Kriechgeschwindigkeit in Baumetallen 
weniger einfach ist als der fiir das Maxwellsche Material vorausgesetzte lineare 
Zusammenhang, kann man erwarten, dass das hier besprochene Modell den Ein- 
fluss von Kriechvorgangen auf die Spannungen in statisch unbestimmten Kon- 
struktionen qualitativ wiedergibt. 


(Received: June 238, 1955.) 


On Stability Questions for Pendulum-Type Equations’) 


By GEorRGE SEIFERT, Lincoln, Nebraska, USA2) 


1. Introduction 


The differential equation 
6+ 40,0) 6= (6), ()= 4) (1 


where / and g have continuous derivatives with respect to 0 everywhere, { is 
continuous in the positive parameter «, g(6) = 0 has simple roots, f(@ + 22, «) 
= {(9,«), and g(@ + 2) = g(0) for all 6, is of interest in certain problems 
relating to the performance of synchronous motors [1] al 

If /(9, %) = f(6), where /(@) > 0 for all 6, conditions have been given 
under which each solution 6(¢) of equation (1) will tend to a finite limit as 
t > +00 [2]; i.e., each solution is stable with respect to a point of equilibrium 

1) This research was supported by the United States Air Force, through the Office of Scientific 


Research of the Air Research and Development Command. 
2) University of Nebraska. 


3) Numbers in brackets refer to References, page 247. 
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of the system described by equation (1). If /(0,«) may have negative values, 
periodic solutions of equation (1) which were excluded by a nonexistence 
criterion of BENDIxson in the previous case-may exist, complicating the 
stability question. It is however possible to extend the notion of stability to 
solutions which coincide with or asymptotically approach some periodic so- 
lution of equation (1); here solutions with period zero (constants) and infinity 
(in the sense to be described below) are allowed. 

Under the conditions on / given in section 2 of this paper we show how a 
value «, may be obtained such that whenever « > «,, each solution of equation 
(1) is stable in the general sense just mentioned. If to these conditions we add 
those on / given in section 3, we may obtain a value «,, such that whenever 
% > a, each solution of this equation is stable in the more restricted sense. 
Under conditions given on / and g in section 4 we show how a value a > 0 
may be obtained such that whenever « < a, this equation will have some 
solutions not stable in either sense. Finally specific examples of functions 
and g are given in section 5 which satisfy all the conditions of the previous 
sections and for which f/ takes on negative values for every « > 0. 

A system equivalent to equation (1) is 


O=y, y=8(6)—/6,a)y. (2) . 


In terms of a (6, y)-coordinate system of a suitably chosen 2-space (which may 
be either the usual cartesian plane or in this case the surface of a circular 
cylinder), the graphs of the solutions of this system are known as its charac- 
teristics or phase trajectories. The singular points of this system are the points 
(6,,, 0) for which g(6,,) = 0; those zeros 0,,; for which g’(@,,;) > 0 are of saddle 
type while those 6,,;,, for which g’(0,;,,) < 0 are either of focal, center, or 
nodal type. We will assume g(0) = 0, g’(0) > 0 which is clearly no restriction 
on generality. 

In terms of the cartesian phase plane a phase trajectory which consists of 
a closed curve corresponds to a periodic solution of finite period of equation (1). 
If a trajectory is closed by the addition of a singular point (necessarily of 
saddle type) it corresponds to a solution which approaches a finite value 6; 
for t + +co and tf > —ov;i.e., may be said to be of infinite period in ¢. 

We now define the two concepts of stability to be used. A solution [6(¢), y(¢)] 
of (2) is said to be G-stable if the positive limiting set (cf. [3], pp. 169-170) of 
its trajectory in the (0, y)-phase plane is either a point or a closed curve. If 
the positive limiting set is a point, the corresponding solution will be said to 
be stable. 

We will refer to the equation 


y'(9) = £9 — 1(6, «) (3) 
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which follows from system (2); the prime here and henceforth denotes dif- 
ferentiation with respect to 0, and 6 will be restricted to the intervalO< 0X22 
unless otherwise stated. 


2. Sufficient Conditions for G-Stability 
We assume that 


a—> +00 


tim (B(0) evel —+00, where (2) = 5— | i(0, «) a0, ‘ 


b(a) 6 — F(8, «) 
gee max sin (6/2) 


|> and F(6, a) =| tH, a) ds. 


Define h,(0,«,¢) = b(«) 0 — F(0,x) + M(a,e) sin (0/2) where M(«,e) = M,(«) + «. 
We first observe that 4,(0, «, €) = ¢ sin(6/2). Next, for « sufficiently large, 
M (a, €) 0 


h’, (0, x, €) + f(8, a) = b(a) “BF SOSts 


M(«, €) | g(8) | g(9) 
eM dies ne - é sin (6/2) = h,(0, a, €) ’ 


i.e., there exists some «,(¢) such that for « > «,(e) we have 


Hee, oc ee 


This last inequality tells us that the slope of any phase trajectory intersecting — 
the graph of y = h,(0, «, €) is less than the slope of this graph at the point of 
intersection for 0 < 6 < 22. This clearly implies that if y(6) is the solution © 
of equation (3) for which 


y(0)=0 and limy’(6)>0, 

6-0, 
then for some c, 0 << c < 2a, y(c) = | 
If we now consider the graph of y =/,(0,«,¢) +a (a> 0), we may simi- — 
larly show that for « > a,(e) the solution y(@) of equation (3) such that 
y(0) = a must satisfy y(0) > (22); we have only to observe that | 


BOWS |g(4) | 
h,(0, a, €) = h,(0, a, €) + a~ 


It follows that for « > «,(e) equation (3) has no non-negative solutions (6) 
for which y(0) < y(2 2). 
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We now define h_(6, a, 2) = b(a) 6 — F(0, «) — Ma, e) sin (0/2), and observe 
that h_(6, «, e) S — e sin(@/2). Since 


hi (0, a, 6) + f(0, a) = bf) — “A ®) ogs 9 


| g(6) | <r | g(9) | ang (9) 
&sin(G/2y == RAO, 6) = A (On, @) ? 


for « sufficiently large, there exists an «_(e) such that for « > a_(e), the slope 
of a phase trajectory intersecting the curve y = h_(0,«,e) is less than the 
slope of the curve at the point of intersection. This implies that if (0) is the 
solution of equation (3) for which 


y(0) =0, limy'() <0, 
then (0) > (22). ; 
Since 


lg(8)| 


h_(6, a, €) 


| g(9) | 


“yO ced) St 


INV 


for a > 0, we have that the graph of the curve y = h_(6, a, e) — a has the 
same properties of contact with the phase trajectories as the curve y = h_(6,«,¢) 
has; this means that for « > «_(e) we can have no non-positive solution (6) 
of equation (3) for which y(0) < (22). 

If we put 

a,(€) = max[a,(e), a_(e)] anddefine a, = inf 0, (€) 
E> 

we have that for « >«, no solution y(0) of equation (3) for which y(0) < (2a) 
is possible. 

The fact that if y(@) is a solution of equation (3) then y(# + 22) is also 
a solution, permits the use as phase space of a right circular cylinder on which 
the §-axis is a circle surrounding the cylinder and having a circumference of 
length 2, and the y-axis is a line on the surface perpendicular to the 6-axis 
at @ = 0. In this phase space the positive limiting set of each trajectory, if 
finite, is one of the following: 

(i) a closed curve surrounding the phase cylinder, 

(ii) a closed curve enclosing a finite region of the cylinder, 
(i) a singular point. 
If « > a, then clearly the limiting set of each trajectory cannot be infinite. 
If it were the first of the above possibilities, it would necessarily be a trajectory 
_ corresponding to a solution y(#) of equation (3) for which y(0) = y(2 7) which 
is impossible. We have then that for « > «, each solution of system (2) is 
G-stable in the sense of the definition given in section 1. 


ZAMP VII/16 
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3. Sufficient Conditions for Stability 


Let 6; (¢ =0,1,..., 2%) be the zeros of g(6) such that 0, — 0, 05, = 2a, 
and 6; ,,. Recalling that g’()) > 0, we have that the 6,, (¢ = 0,1,..., ) define : 
the saddle singularities of the system (2). We assume that for each 7 = 0, 1, 
2,..., 2 — 1 there exists an interval . 


ay; < 606, such that! 03, S a,7 ss Ofpey ee, Saya 


in which 
G) ((0;.«) 230 “for “age 0 = ey 


% 


(ti) Ags (a) MY (a) >+00 as a>-+00, 


Doi42 — ei 


(iit) By ,(«) — = (4) ++ co as a ++ co; 


a =a 
here, 
995 +2 


: th /(0, 2) 40, 


Aai(a!) os Oei42 — a2; 
Oy 


| 
(0 — ag;) Ag;(«) fin u, «) du z 
M§) (a) = max 924 
99,5950; ,5 3 0 — ay; , 
[sn oa] | 
bo4 
Bs (a) = — a | HKG, &) as 
boi 69; J 
6 
(0 — 63;) By i(a) — [iu, a) du 


bo 
S 
+ 


92550 bo; 


My}() =, min 


We will restrict ourselves to the case i= 0 and point out that for the 


other values of 7, a procedure identical, except for subscripts, with this one 
will apply. 


Define the function 


6 
HS(0, a, &) = (0 — ag) Ag(a) — f f(u, a) du — [M$ (a) + e] sin [pata 


Vol. VII, 1956 On Stability Questions for Pendulum-Type Equations 243 


We observe that for a < @ < 64, 


WY) (0, «, €) S—esin een } 
0 


that hj? (B2, %, €) = hi (a, «, €) = 0, and that in view of condition (ii) given 
above, there exists an «_(e) such that for « > «_(e) and a <0 < 4, 


h'(0, «, €) > BA Baal 6, a) . 
0 ( ) nie) (0, w: e) I( , a) (6) 
If a) > 6) we have, since g(a) > 0, this inequality at 6 = a, in the sense 
that 
ies >—oco as $>a,;. 
AW) (0, x, €) 


Inequality (6) means that each trajectory intersecting the graph of 
y= he(0, ass) ih “ays 0 3.0, 


has a slope less than the slope of the graph at the point of intersection. Hence 
the solution y,(@) of equation (3) for which 


* 


yo(0s) =0 and lim y6(8) > 0 


6—6; 


‘must have a zero in the interval a, << 6 < 6,;1.e., the trajectory corresponding 
to this solution must, after emerging from the point (6, 0) with increasing 7, 
intersect the 6-axis in the interval ay < 0 < 6. Let 4 be the smallest value 
of ¢ for which y,[6(¢,)] = 0. We will now show that a < 0(f) S 4). 

For suppose 6, < 6(t) < 6,. We have from equation (3) that the trajectory 
corresponding to this solution can have no tangents parallel to the y-axis in 
this interval. By the mean value theorem for derivatives, one can determine 
a monotone decreasing sequence {&,}, E, > O(f), such that 3(é,;) < 0. Since 
e(E;) > g[6(%)] < 0, and y.(&,) > 0, while y,(&;) < 0, there exists an 1 suf- 
ficiently large for which 


, =~ B(En) = 
Vo(En) = Vo(En) flEgra) 20% 


a contradiction. Hence we must have a < O(f) S 6,. We now define the 
function 


6 
HO, aw, €) = (0 — 9) Bola) — /[ f(t, «) du —[M!(«) — e] sin| — 2 x] 
% 
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and observe that hi”) (Oy, «, €) = At” (bp, «, €) = 0, that 


1(0, a, €) = esin[ cata for 0<0< by, 
(iy 0 


and that in view of condition (iii) given above, there exists an «,() such that . 


for a > «,.(€) 


h® (0, a, €) > — ML {(0,«) for 650< by. (7) 
0 


for 0 > 0, . 
Inequality (7) means that each trajectory intersecting the graph of 


=h(0,0,¢6) in 0505 by 


has a slope less than the slope of the graph at the point of intersection. Hence ~ 


the trajectory emerging from the point (8, 0) with positive slope and increasing 
¢ intersects the 6-axis first in the interval 0, < 6 < bg; as in the previous case, 
the interval 0) < 0 < 6, can be excluded. 


We now set «(e) = max [«_(e), «,(e)] and consider « > «'®)(e). Clearly, any 


closed trajectory contained in the strip 0) < 0 < 6, must necessarily be con- 


tained in the strip a) << 6 < bo. Since /(8, a) > O in a < 0 < dy we have by ~ 
the non-existence criterion of BENDIxsoN that no closed trajectory in this — 


strip, and hence, in 0) < 0 < 6g, is possible. We also observe that no trajectory 
which may become a closed curve upon addition of either (05,0) or (6, 0) 
is possible. 

We now proceed similarly to obtain «(s) (¢ = 1, 2,...,—1) such that 
for a > a'(e), no closed trajectory is possible in the iaicual 03 S054 
Let us now assume that a closed trajectory exists in the phase plane for 


ot 2 Og(é) = max {«'(e) (e)}. Clearly, this implies the existence of a closed tra- | 
jectory enclosing some saddle point (0;, 0) (Rk <1) and intersecting the 0-axis — 


at 8, < 95; such that for some; < k, 02; < 0, < 05; , 1. However, for a > as(e), 
this trajectory intersects the graph of y = hf ’}(0, «, €) with a greater slope than 
the slope of the graph. This is impossible, and we conclude that for a > a,(€), 


no closed trajectory is possible anywhere in the phase plane. 
If we put 


As == inf %»(€) 
e>0 


we have that for « > a, any finite positive limiting set of a solution of (2) 
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must be a point; i.e., any solution of (2), if G-stable, must be stable. Finally, 


if x > max(«,, %), where the «, is defined in the’previous section, each solution 
of (2) is stable. 


4. Conditions for the Existence of Nonstable Solutions 


We say that a solution of (2) is nonstable if it is not at least G-stable, and 
assume here that 


five a)|d9>0 as a>0, (i) 
[eo do>0. (ii 


0 


Condition (ii) above assures the existence of at least one zero 6, of g(6) such 


that g’(4)) > 0 and : 


if g(u) du > 0 
Oy 
for 0,< 65 0,+ 22 [4]. 
Consider the solution y9(@) of equation (3) for which 


Yo(9o) =O and dim, 96(4) = Org 


Suppose that y,(#) has zeros in 6) < 6 < 6)+ 22 and denote by , the least 
such zero. As was observed in a previous section, we have 0, S< << 0+ 2a, 
where 6, is the smallest zero of g greater than 6). From equation (3) we have 
that y,(9) must satisfy 


¥o(9) yo(8) = g(8) — £(6, a) vo(4) - (8) 


Integrating equation (8) from 6, to 6, we have 


9 4 


[s(0) 40 = | 100, 2) v(6) a0. (9) 
9 


Let 
max ¥o(8) = Yo(9m) = Ym - 


6,<6<4 


Integrating equation (8) from 4, to 6,, we have 


Om Om 
Yn, = | (6) do — | {(9, ) yo(9) 40 < Gy + F(@) Ym» 


0 % 
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where 


27 


Ge= ) max [e) ds, and F(a) = | |H6, a) | ds. 


9 
Oy 8S eg: a 


Hence ym < F(«) + [F2(«) + 2 Gi]? = Ko(«), which used with (9) gives 


[8(0) 40 < Kola) Fle) (10) 


Since 


; 

there is m, > 0 for which 
3, . 

[8(0) 40 = mo. . 

6, : 

Since as « >0, F(x) >0 and K,(«) > (2 G)'* we may determine an «(O ) © 


such that whenever 0 < « < a(Oy) we have 
% 
- [e(0) a0 = my > Kola) Fle) 
8 ¢ 
which contradicts (10). 
We thus conclude that for « < «(0p), Y9(9) can have no zeros in 


6,<0< 6,122; i.c., 9 (0,+22) =O. 


We now refer to cylindrical phase space and observe that the trajectory for — 
any solution [0(¢), y(¢)] of (2) containing (5, Yo) where yy) > 0 can have a 
closed curve as positive limiting set only if the closed curve surrounds the 
phase cylinder. Hence, such a solution cannot be G-stable and we have shown - 
that for a < «(05) such een as exist. 
If we now put % = max{«(0,)} where {6; pi is the set of zeros of g(6) in 
0 <0 S22 for which 
D 
[e0)d>0, 0>6,, 


6; 


and «(6;) is determined as was a(Q)) above, we have that for « < a, solutions 
of (2) exist which are not stable. 
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5. A Special Case 
We consider the case where 


ea) 2 i 
g(6) = sin—[—* + sinZ, (0,0) =a (1-+ bsin2 6) . 


We observe that 6)= 0, 6,= 32/2, 0,=2m. Using the terminology of 
section 2 we have that b(«) = «, 
| bsin®@ ad 
Myla) =, max |e = 83-5. 
and hence b(x) — Mg(«)/2 = «(1 — 4 V3 b/9) > +00 whenever b < 3 3/4. 

Taking a) = a and by) = 13 2/8 the conditions given in section 3 are satis- 
fied for some 6 > 1; we omit the details. 

The conditions of section 3 are clearly satisfied for any b. 

We point out in conclusion that apart from a translation of magnitude 
52/4 in @ the /(0, «) of this special case is precisely of the same form as the 
one considered by EDGERTON and FouRMARIER in the treatment of a problem 
concerning the performance of salient-pole synchronous motors [1]. 
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Zusammenfassung 


Die Lésungen 4(¢) der Differentialgleichung 


d?6 do 

pet a 

werden untersucht, worin /(0 + 22, «) = f(0, «), g(0 + 2 2) = g(@), g(@) einfache 
Nullstellen besitzt und /(6, «) fiir positive « negative Werte annehmen kann. 
Bedingungen werden angegeben, fiir welche jede Losung beschrankt bleibt, wenn 
t-> + co und sich asymptotisch einer periodischen Lésung annahert; diese pe- 
riodische Lésung kann auch eine Konstante sein oder die Periode + oo haben. 
Bei engeren Bedingungen gilt fiir jede Losung 6/(?), dass 


_ jim 6 (4) = (hie 


wo 0, eine endliche Konstante ist. Schliesslich wird ein Beispiel aus der Elektro- 
technik angefiihrt, welches allen diesen Bedingungen genigt. 


(Received: June 11, 1955.) 
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Kurze Mitteilungen — Brief Reports - Communications bréves 


Biegeschwingungen verwundener, einseitig eingespannter 
und am andern Ende gelenkig gelagerter Stabe 


Von Max ANLIKER, Brooklyn, New York, N. Y., USA?) 
Unter dem Begriff « Verwundene Stabe» versteht man Stabe, die im unge- 


spannten Zustand wie ein Propellerblatt oder eine Turbinenschaufel verdreht 
sind (vgl. Fig. 1). 


Figur 1 Figur 2 
Sttick eines spannungslos verwundenen Einseitig eingespannter und am andern Ende 
Stabes. gelenkig gelagerter Stab. 


_Die Biegeschwingungen (Standschwingungen) von solchen Staben werden 
bei der blichen Beschrankung auf von erster Ordnung kleine Gréssen durch die 
Torsionsschwingungen nicht beeinflusst und kénnen daher separat behandelt 
werden. Die Eigenkreisfrequenzen x lassen sich mit gewissen Einschrankungen 


L D . . , 3 
ae shee 2 of Aeronautical Engineering and Applied Mechanics, Polytechnic Institute of 


3 f 
) A. Tr6scu, M. ANLIKER und H. Zrecurr, Quart. appl. Math. 12, 163 (1954). 
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exakt berechnen, wie A. Tr6scH, M. ANLIKER und H. ZIEGLER?) in einer kiirzlich 
erschienenen Arbeit gezeigt haben. Lasst man namlich nur homogene isotrope Sta- 
be zu, die einen konstanten Verdrehungswinkel pro Langeneinheit (w) aufweisen, 
und nimmt man zudem an, dass die Masse pro Langeneinheit (uw), sowie die 
Biegesteifigkeiten (x, 8) beziiglich der Hauptachsen konstant sind, so ergibt sich 
fiir die Bewegung der elastischen Linie-ein System von 2 linearen, homogenen, 
partiellen Differentialgleichungen vierter Ordnung. Mit dem Separationsansatz 
und 8 linearen homogenen Randbedingungen fiihrt die Lésung der Schwingungs- 
aufgabe auf ein Eigenwertproblem achter Ordnung, dessen Eigenwerte mit Hilfe 
eines modernen Rechenautomaten berechnet werden kénnen. 


: | 
65 1 SN 


ON 


v= 1/256 
50 ma pore 
WW 
45 \ L. 
\ O 
40 eae 
\ y=1/25| 7 ae 
35 aos eet 
30 | , 
v=1/4 
25) = hye a 
at 


15 : 
0 30 60 90 120 150 180 OT ng30" "G8 

Figur 3 Figur 4 

Die ersten 2 Frequenzkurven m?(®) fiir 3 verschiedene Verhaltnisse der Biegesteifigkeiten. 

Bei rechteckigem Querschnitt entsprechen den gegebenen y-Werten die Seitenverhaltnisse 1:2, 


1:5 und 1:16. Kreisfrequenz x; = m3 /P Vo) hia 
(J Stablange, « Biegesteifigkeit, Masse pro Langeneinheit, ® totaler Verdrehungswinkel). 


20 
° 
VET eres 
—<— 
a’ 90" 1205 150 4718 


(0) 
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In der oben erwahnten Publikation wurden so die ersten 4 Eigenkreisfrequen-_ 
zen fiir den einseitig eingespannten Stab mit ideal schmalem Querschnitt?) 
exakt ermittelt. Dazu wurde die programmgesteuerte Rechenmaschine des Insti- 
tutes fiir angewandte Mathematik der ETH. verwendet. Im Anschluss an diese 
Arbeit werden hier die Resultate angegeben, welche mit Hilfe des gleichen Rechen- . 
automaten beim einseitig eingespannten und am andern Ende gelenkig gelagerten 
Stab (vgl. Figur 2) bei verschiedenen VerhAltnissen v = «/f der Biegesteifigkeiten — 
gefunden wurden. 

In den Figuren 3 bis 8 ist der Verlauf der Frequenzgrésse m? bzw. 


ee ; 
m= (1 Lange des Stabes) 


— 


, 
4 
: 
Figur 5 ; 
Dritte Frequenzkurve m3(®) fiir 3 verschiedene v-Werte. 


) Das heisst eine der Bie ifigkei i 
Sfcsnl aaa gesteifigkeiten wird als unendlich gross, die andere als endlich 


Vol. VII, 1956 Kurze Mitteilungen — Brief Reports — Communications bréves 251 


fiir 3 verschiedene Werte von y als: Funktion des totalen Verdrehungswinkels 
gy aufgetragen. Dabei sind durch kleine Kreise die exakt berechneten Punkte 
gekennzeichnet. Die gestrichelten und strichpunktierten Linien stellen Naherungs- 
lé6sungen dar und wurden rein analytisch mit Hilfe von Reihenentwicklungen 
ermittelt. 

Wie aus diesen Figuren hervorgeht, hangt der Verlauf der Frequenzfunktionen 
m(qg) insbesondere im praktisch interessanten Bereich 0° < pm < 45° stark von 
der Verwindung ab. Ausserdem zeigt sich der Einfluss der Schlankheit des Quer- 
schnittes recht deutlich. So steigt zum Beispiel m,(q) bei kleinen g um so steiler 
an, je schmaler der Querschnitt ist, wahrend m,(y) mit abnehmendem y zuerst 
ansteigt, dann abfallt. Bei grossen stellt man fest, dass je zwei aufeinander- 
folgende Frequenzkurven m,,;_, und m,; der gleichen Asymptote zustreben‘). 

Der Verfasser ist zur Zeit mit der Untersuchung weiterer Falle (das heisst 
anders gelagerter Stabe) beschaftigt und beabsichtigt, die zugehdrigen Resultate 
spater zu verdffentlichen. 

Bei dieser Gelegenheit méchte der Verfasser dem Direktor des Institutes fiir 
angewandte Mathematik der ETH, Herrn Prof. Dr, E. STIEFEL, sowie dem Pra- 
sidenten des Schweizerischen Schulrates, Herrn Prof. Dr. H. PALLMANN, fiir die 
wohlwollende Unterstiitzung und Férderung dieser Arbeit bestens danken. 
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Figur 6 
Die ersten 4 Frequenzkurven fiir vy = 1 /4. 


4) Hine ausfiihrliche Darstellung dieser Resultate findet sich in der von der Eidgendssischen 


Technischen Hochschule genehmigten Promotionsarbeit Biegeschwingungen verwundener, einseitig 
eingespannter und am andern Ende gelenkig gelagerter Stéibe, vorgelegt von M. ANLIKER. 
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0 90 180 -270 360 450 540 630 720 8710 
Figur 8. Die ersten 3 Frequenzkurven fiir » = 1 /256. 
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Summary 


: The lateral and torsional vibrations of twisted rods can be treated separately 
if we consider as usual only first order terms. The eigen-frequencies of the lateral 
vibrations can be calculated exactly if we restrict ourselves to isotropic homo- 
geneous rods with constant mass and twist per unit length and constant principal 
flexural rigidities. In this paper the eigen-frequencies for a rod built in at one 
end and supported at the other are given for 3 different cross-sections. 


(Eingegangen: 9. November 1955.) 


Der symmetrische Kardankreisel unter einem Moment an der Achse 
des dusseren Rings 


Von Hans ZIEGLER, Ziirich!) 


1. Problemsteliung 


Einer der bekanntesten Kreiselversuche besteht darin, dass man versucht, 
den ausseren Ring eines sonst kraftefreien, kardanisch gelagerten Kreisels um 
seine Achse zu drehen. Solange die Figurenachse des Kreisels nicht mit der Dreh- 
achse dieses Rings zusammenfallt, kommt die angestrebte Bewegung nicht zu- 
stande; dagegen dreht sich der innere Ring im dusseren, bis die Vektoren, welche 
die Winkelgeschwindigkeiten der Eigenrotation und der angestrebten Drehung 
darstellen, gleichgerichtet sind. 

Dieser Effekt, der sich mit Hilfe des Drallsatzes qualitativ leicht erklaren 
lasst, scheint bisher quantitativ nicht untersucht worden zu sein. Er lasst sich 
aber unter der Annahme, dass der Kreisel schnell sei, elementar behandeln, vor 
allem, wenn die Massen der Kardanringe, ihre Elastizitat und die Lagerreibung 
vernachlassigt werden. 


2. Bewegungsdifferentialgleichungen 


Unsere Abbildung zeigt schematisch einen symmetrischen Kardankreisel mit 
den Eulerschen Winkelgeschwindigkeiten p= w,, 8 und »= , unter dem 
Einfluss des an der Achse des dusseren Rings angreifenden Momentes M. Das 
raumfeste Koordinatensystem (%, y, 2) ist so gewahlt, dass die z-Achse, um die 
sich der aussere Ring dreht, vertikal nach oben weist. Die ¢-Achse des korper- 
festen Systems (&, 7, ¢) ist als Figurenachse so orientiert, dass sich der Rotor 
im Rechtsschraubensinn um sie dreht, mithin 


wo, > 0 (2219 


ist. Die Drehung des inneren Rings relativ zum dusseren erfolgt um die Knoten- 
achse x; der Richtungssinn derselben wird beim Kardankreisel zweckmassig ein 
fiir allemal (wenn auch willkiirlich) fixiert, so dass dem Winkel @, der mit x zu- 
sammen der Rechtsschraubenregel geniigt, das Intervall —7 = # = a zugestan- 
‘den werden muss. Die Knotenachse x kann mit der Querachse 4 und der Figuren- 
achse ¢ zu einem dritten Bezugssystem (x, A, ¢) zusammengefasst werden. 


1) ETH. 
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Da der Schwerpunkt des Kreisels mit seinem Drehpunkt zusammmenfallt und 
die Ringe als masselos vorausgesetzt werden, setzt sich das resultierende Moment 
am Rotor aus dem Moment M mit vertikalem Vektor und einem Moment M” mit 
horizontalem Vektor zusammen, das von den Reaktionen in den Lagern des } 
ausseren Rings herriihrt. Zerlegt man es im System (x, A, €), dann miissen die 
Komponenten M, und Mz verschwinden, da die Lager auf der w- und ¢-Achse 
als reibungsfrei vorausgesetzt sind. Die Momente M und M’ miissen sich also . 
zum resultierenden Moment M) in der Querachse zusammensetzen lassen. Hieraus 


Momente am Kardankreisel. 


folgt, dass M’ in der Querebene (A, €) legt und 


M M 


ie Mame Te contig Rode (os Ga. 
) 


ist; die Winkel #= 0 und #@= + aw sind dabei auszuschliessen. 
Wird der Kreisel als schnell, also 


| op |< op, |b1<e, (2.3 


§ 
‘ 
| 
angenommen, dann kann der Drallvektor D mit Hilfe des axialen : 


momentes C und des Einheitsvektors ¢in der Figurenachse naherungsweise durch 
D=Co,é (2.4) 

dargestellt werden. Nach dem Drallsatz ist dann 
Co,é@é=M,, Co,=Me, (2.5) 


wobei M, die Projektion des resultierenden Momentes auf die (x, 4)-Ebene, das 
heisst den Vektor mit den Komponenten M,,, M, darstellt. Der Vektor z, der 
als Geschwindigkeit der Kreiselspitze gedeutet werden kann, hat (Abbildung) im 
System (x, 4, ¢) die Komponenten w, sin #, — #, 0. Man kann also statt (2.5) auch 

C © Sind = My, —Co,6=Mi, Co,= Mz (2.6) 


schreiben und erhalt hieraus, da M, sowie M ¢ Null sind und M, durch (2. 2) 
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gegeben ist, 
M 


@, == Oi; Gereee eee e= 
p Ca,sin® Oo, const . (2.7) 


3. Integration 


Nach (2.7) kommt beim schnellen Kreisel unter dem Moment M keine Priizes- 
sion zustande; die Eigenrotation bleibt gleichformig, und die Figurenachse hebt 
oder senkt sich, je nachdem das Moment M positiv oder negativ ist, das heisst 
eine Rechts- oder Linksprazession anstrebt. 


Die zweite Beziehung (2.7) lasst sich ohne weiteres integrieren. Fiihrt man 
namlich mit 


cos} = w (3.1) 


die z-Koordinate der Kreiselspitze ein, dann geht sie in 


w= Biches Se 
One bag) 
liber; es ist also 
t 
1 H 
w= Wy + Co, | Me. | (3.3) 
0 


wenn mit w,) = w(t = 0) die anfangliche z- Koordinate der Kreiselspitze bezeichnet 
wird. 
Fail a: Setzt man zum Beispiel 
M = M, | sin? | (3. 4) 


an, dann folgt aus der zweiten Beziehung (2.7) 
aie Ne er, <a (B20): (3. 5) 


man erhalt also insbesondere fiir die volle Schwenkung der Figurenachse die 
Dauer 


Cw 
w= 7 re. (3.6) 
|M,| 
Fall b: Nimmt man 
Diem VA (3.7) 
konstant an, dann liefert (3. 3) 
M : 
bet ae Seer ((3. 8) 
und damit insbesondere die Schwenkdauer 
ysuiacpeceeer '(3. 9) 


|Mo| 


Den zweiten Beziehungen (2.7) und (2.3) zufolge kann freilich der Kreisel mit 
nahezu vertikaler Figurenachse nur dann als schnell gelten, wenn — wie 1m Falle 
a—|M| mindestens wie | sind? | gegen Null geht. Trifft dies - wie im Falle b — 
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nicht zu, dann gilt (3. 3) nur in einem reduzierten Winkelbereich « <|0| Sa —€, 

und fiir| #|< esowie| #| > a — emiisste die Bewegung mit der strengen, auch 

fiir nicht-schnelle Kreisel giiltigen Theorie ermittelt werden. Falls aber der 
Kreisel bei annihernd horizontaler Figurenachse geniigend schnell ist, fallt e sehr 

klein aus, und die im Falle b aus (3.9) folgende Zeit #* kann damn als gute Nahe- 

rung fiir die Dauer einer vollen Schwenkung, das heisst fiir praktische Zwecke 
doch wieder als volle Schwenkdauer angesprochen werden. 


Summary 


The motion of a high speed gyroscope under the influence of a couple applied 
to the external Cardan ring is discussed. 


(Eingegangen: 16. Februar 1956.) 


Measurement of the Thermoelectric Power of Lead Sulphide 


By Art Art AraFa, IBRAHIM SAVED SHAFIE, and FaTHI SULTAN AHMED SULTAN, 
Alexandria, Egypt?) 


Abstract 


The thermoelectric power of natural specimens of lead sulphide crystals has 
been measured between room temperature and 700°K. Specimens which are of 
the n-type remain as such throughout the temperature range. The p-type spe-— 
cimens, however, show an initial small rise followed by a gradual decrease and a 
subsequent change of sign at a temperature which, in some specimens, is well 
defined. Further increase of temperature results in a negative extremum. From 
this stage onwards there seems to develop a common behaviour for both types 
of specimens. Some crystals which were initially of the p-type changed over to 
the n-type and remained as such throughout subsequent investigation. Among 
the crystals one showed a metallic conduction and an extraordinary large thermo- 
electric power. 


Introduction 


A good deal of literature is now available about the thermoelectric and other 
electrical properties of semiconductors. Recently much work has been devoted to 
the thermoelectric properties of germanium and silicon and a number of metal 
sulphides, selenides and tellurides. LarK-Horovitz et al.[2]?) obtained the 
thermoelectric power curves of aluminium-doped germanium samples with dif- 
ferent carrier densities. MippLeton and ScaNLon [3] investigated aluminium- | 
doped (p-type) and antimony-doped (n-type) germanium samples. The sign of 
the thermoelectric power in the low temperature or impurity range is the sign - 
of the carrier, i. e., the p-type samples have positive thermoelectric power and 
the n-type negative power. The thermoelectric power of an n-type sample 
remains negative throughout. Beginning at the low temperature end of the 
curve, the thermoelectric power increases in magnitude with rising temperature, | 

q 


1) Physics Department, Faculty of Science, Alexandria University. 
2) Numbers in brackets refer to References, page 264. 
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passes through an extremum and then decreases in magnitude, approaching a 
curve common to all samples. The behaviour at high temperatures is so similar 
that it may be considered as characteristic of the intrinsic material, Beginning 
at the low temperature end for a p-type sample, the positive thermoelectric 
power rises to a maximum, then decreases fairly rapidly, eventually changing 
sign. This is followed by an extremum-value and a gradual approach to the 
intrinsic range. 


PUTLEY [4] carried out similar work on lead selenide and telluride and obtained 
qualitatively the same behaviour. 

The present work deals with experimental technique and results of measure- 
ments of the thermoelectric power of natural crystals of lead sulphide and its 
variation between room temperature and 700°K., 


Experimental Procedure 


The crystals were cleaved into small rectangular blocks of approximate size 
4x4x2mm. The specimen is enclosed between two nickel sheets which at one 
end are connected to two copper-constantan thermocouples (Figure 1). At the 
other end they are connected to a multi-range Habicht galvanometer M. Pressed 
against the specimen and mica-insulated therefrom are two small hot plates H 
and H’ which are electrically heated by two small variable A. C. voltages. A 
two-way key A B enables in position A a direct connection to a Tinsley potentio- 
meter Type 3387 B which reads down to 1 wV. Position B enables a measurement 
of the resistance of the specimen. Thus a simultaneous determination of thermo- 
electric power and conductance is possible at different temperatures of the crystal. 
The temperature of the crystal is taken as the mean of the temperatures of the 
two faces, the difference being kept small (12°). The thermoelectric power is 
obtained by dividing the potentiometer reading d6 by the corresponding difference 
of temperature dT. 


Figure 1 


Electrical setup for measuring thermoelectric power. 


Results 


Some specimens begin, in the low temperature or impurity range, as p-type, 
the carriers being predominantly positive holes. Figure 2 illustrates the behaviour 
of this type of specimen. The positive thermoelectric power d0/dT shows an 
initial small rise to a maximum of 0-27 to 0-35 mV/degree at about 380° K, followed 


ZAMP VII/17_ 
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by a steady decrease. This indicates the onset of mixed conduction. As the tem- 


perature increases electrons become progressively the more dominant carriers and 


a change over in the sign of the thermo-electric power takes place. In three con-— 


secutive runs (1), (2) and (3) the thermoelectric power is measured at temperatures 


ranging between room temperature and about 700° K. The crystal is left to cool — 
down to room temperature after each run. It is noted that the change of sign takes © 


‘ 


place at almost the same temperature of about 550°K. As the temperature is — 


further increased d0/dT reaches an extremum at about 650°K, after which it 
decreases in magnitude. 


0-4 
03 


0-2 


dO/d/T in mV/degree 
° 
° ms 


\ 
2 
EN 


-02 


Figure 2 


Crystal S-6. Temperature variation of the thermoelectric power of a p-type specimen. (1), (2) and (3) 
are three consecutive runs. After each run the specimen is left to cool down to room temperature. 


Figure 3 illustrates the behaviour of an n-type specimen S—7. The thermo- 
electric power remains negative throughout; indicating a dominance of the 
electrons as carriers. Four consecutive runs are shown. In the first run (1), d0/dT 
decreases slowly in magnitude with rising T. In the following runs (2), (3) and (4), 
however, d0/dT increases in magnitude rather rapidly to an extremum at a 
temperature lying between 425° and 470°K. This indicates again the onset of 
mixed conduction due to thermal release of electrons from and holes in the full 
band. Further increase of temperature leads to the intrinsic range. The fact that 
d0/d7 never changes sign in this type of specimen is due to the fact that the 
number of electrons is, at all temperatures, greater than the number of the less 
mobile holes. 

Figure 4 refers to the same n-type specimen S—7. It shows the relation between 
the logarithm of conduction o and the reciprocal of the absolute temperature. 


ee a 


ele initial thermoelectric and conduction behaviour of the specimen, as indicated — 
y curves (1) in Figures 3 and 4 respectively, seem to obey the laws of classical _ 
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Figure 3 


Crystal S-7. Temperature variation of the thermoelectric power of an n-type specimen. (1), (2), 
(8), and (4) are three consecutive runs. 


semiconductor theory. The activation energy ¢« as calculated from the slope of 
the logo versus 1/T line (1) is about 0-156 eV. It follows that the chemical 


potential, or free energy per electron, in this specimen is 
éE 
= ae Bae 0-078 eV . 


Substituting this value in the classical expression for the thermoelectric power 
of an extrinsic semiconductor: 

Oe ok ( pase: e 

dT e kT 


where & is BoLttTzMANN’s constant (0-68 x 10-4 eV/degree), we have at room 
temperature 

te = — 0:86 10-4 (2 + pore) V/degree = — 0-44 mV/degree 
which is in good agreement with the observed room temperature value —0-47mV 
per degree [Figure 3 (1)]. 
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Figure 4 
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Crystal S-7. Logarithm of conductance against 1/T in the same consecutive runs of Figure 3. The ~ 


activation energy is about 0:16 eV. In the following runs it is almost negligible. 


The rather abnormal behaviour exhibited by the nonlinear relations between 
the logarithm of conduction and the reciprocal of the absolute temperature 
[Figure 4 (2), (3), (3) and Figure 5 (1), (2), (3), (4), (5)] cannot be explained on the 
basis of the simple Wilson-Gudden model of a semiconductor, according to which 
such a relation is linear. The minimum of conductance [Figure 4 (4)] may be 
described as the boundary between two ranges: (a) the low temperature range 
where the conduction is determined almost entirely by the high mobility of the 
relatively small number of carriers, and (b) the high temperature range where 
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Figure 5 
Crystal S-4, Logarithm of conductance against 1/T in five consecutive runs. 


the mobility is much reduced and the conduction is determined by the large 
number of carriers. The maximum of conductance (Figure 7) is, however, difficult 
to explain on a similar basis. The appearance of maxima or minima or both in 
the conductance-temperature graphical relations as well as the ‘metallic’ beha- 
viour of certain semiconductors seem to find explanation in a semiconductor 
model which was first proposed by Giso.r [1]. He assumes that deviations from 
stoichiometric composition, which are pronounced especially in the oxides, sul- 
phides, selenides and tellurides of the multi-valent metals, give rise to donor levels 
whose energy lies above the bottom of the conduction band, or to acceptor levels 
whose energy lies below the top of the highest full band. The energy values asso- 
ciated with these impurity levels are unimportant, since the impurity atoms 
giving rise to them are almost completely ionized at all temperatures. Just below 
the conduction band GisoLF assumes a continuous band of levels into which free 
electrons may be trapped. The material is a typical semiconductor, i. e., possesses 
a negative temperature-resistance coefficient if the density of impurity donors 
(already ionized) is so much less than the density of electron traps that the 


262 Kurze Mitteilungen — Brief Reports - Communications bréves ZAMP 


effect of temperature is a net increase of the number of free carriers contributed — 


either from these traps or from any impurity levels that happen to lie below the 
conduction band. This condition is satisfied by most semiconducting materials. 
On the other hand if the density of the ‘Gisolf impurities’ is dominantly large, 
then the semiconductor is metallic. In this case the effect of temperature on the 


concentration of free carriers is neglected in comparison with its effect on their — 


mobility. Under these conditions the impurities are so close together that appre- 
ciable overlapping of the impurity wave functions occurs and the impurity atoms 
play much the same role as the electrons and ions in a metal, and the excess elec- 
trons form a degenerate electron gas which moves through the irregularly distri- 


buted array of positive ions. As the intrinsic range is approached, however, the ~ 


300 400 500 600 7K 


AG/d'T in mV/degree 


Figure 6 


—5. Temperature variation of the thermoelectric power of a metallic specimen. d0/dT is 
large and increases rapidly with T. 
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Crystal S—5. Logarithm of conductance against 1/T for the metallic specimen. 
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Figure 8 
Crystal S-8. Irreversible change from p- to n-type after the first run. 
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number of the carriers and the conductivity begin to increase with increasing 
temperature. : 

There is evidence for the presence of impurity levels lying either just below 
or in the conduction band. It follows that the energies required to remove 
electrons from such impurities are so small that a negligible fraction of the excess 
electrons are bound to their impurities at room temperature. Refering to Figure 4, 
the initial room temperature behaviour as indicated by the curve S—7 (2) is almost 
temperature independent. This may be explained as due to a negligible activation 
energy; the small increase in the number of carriers due to increase of temperature 
being counterbalanced by the decrease in mobility. 

Figures 6 and 7 illustrate the thermoelectric behaviour and conduction 
respectively with temperature of the specimen S—5. This specimen shows a metallic 
behaviour and, moreover, possesses a very high thermoelectric power lying be- 
tween 2 and 3 mV/degree at a temperature of about 650° K. 

The irreversible change over from p- to n-type exhibited by specimen S—8 
(Figure 8) is due to the romaval of acceptor impurities, e. g., evaporation of 
excess sulpher. Heat treatment to about 700° K causes this type of impurity to 
disappear. 

The complete analysis of these results and quantitative determination of 
mobility and concentration of carriers at different temperatures will be possible 
when further experimental work on the Hall constant and conductivity of the 
same material is carried out. This work is in progress. 
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P Zusammenfassung 

_ Die thermoelektrische Kraft von natiirlichen Proben von Bleisulfid wurde von 
Zimmertemperatur bis 700°K gemessen. Proben vom n-Typ blieben iiber den 
ganzen Temperaturbereich unverandert. Proben vom -Typ zeigten jedoch eine 
Anderung des Vorzeichens. Eine weitere Erhdhung der Temperatur fiihrte zu 
einem negativen Extremwert der thermoelektrischen Kraft. Von dieser Tem- 
peratur an aufwarts scheint sich ein gleichartiges Verhalten der Proben beider 
Typen zu entwickeln. Einzelne Kristalle, die anfangs dem p-Typ angehorten, 
wechselten in den x-Typ iiber und behielten diese Eigenschaft wahrend der ganzen 
folgenden Untersuchung bei. Von den Kristallen zeigte einer eine metallische Leit- 
fahigkeit und eine ausserordentlich grosse thermoelektrische Kraft. 


(Received: January 12, 1956.) 
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Conference on Nuclear Structure 


The University of Pittsburgh announces its third Conference on Nuclear Struc- 
ture, to be held in Pittsburgh during June 6, 7, and 8 of 1957. Communications 
regarding the conference should be addressed to N. AUSTERN, at the Sarah Mellon 
Scaife Radiation Laboratory, University of Pittsburgh, Pittsburgh 13, Pennsyl- 
vania, USA. It is expected that this conference will be held regularly thereafter 
at two-year intervals. N. AUSTERN 


Jahrestagung GAMM 1957 


Die nachste wissenschaftliche Jahrestagung der GAMM (Gesellschaft fiir Ange- 
wandte Mathematik und Mechanik) findet vom Dienstag, dem 23., bis Samstag, 
dem 27. April 1957 an der Universitat Hamburg statt. Die drtliche Tagungslei- 
tung liegt in Handen von Prof. Dr. Kart WieGHarp (Universitat Hamburg). 

J. HeEINHOLD 


Buchbesprechungen — Book Reviews — Notices bibliographiques 


Theorie der linearen Wechselstromschaltungen. Von WILHELM CAUER. 
2. revidierte Auflage, herausgegeben und aus dem Nachlass erganzt von WILHELM 
KLEIN und Franz M. Petz (Akademie-Verlag, Berlin 1954). 769 S., 461 Abb.; 
DM 48.-. 

Es ist sehr erfreulich, dass das Werk des 1945 kurz vor dem Ende des Krieges 
ir. Berlin tragisch ums Leben gekommenen Ingenieur-Mathematikers eine Neu- 
auflage erfahrt. Bei den Ingenieuren galt die erste Auflage als schwer lesbar. Sie 
enthielt allerdings Ideen und Entwicklungen, fiir welche die Technik erst heute 

-reif ist. Deshalb wird die Neuauflage durchaus modern wirken, ja erst wirklich 
fruchtbar sein. Die Herausgeber haben es iibrigens verstanden, bei pietatvoller 
Respektierung des Buches von CaveEr, durch Umstellungen in der Kapitelreihen- 
folge, Einfiigen von Untertiteln, Kleindruck von nicht absolut notwendigen 
Beweisen, erlauternden Fussnoten und Literaturhinweisen bis 1953, den Inhalt 
dem Leser naherzubringen. In der Netzwerktheorie unterscheidet man Analyse 
und Synthese. In beiden Richtungen wirkte CavER bahnbrechend. Die Revision 
bezog sich hauptsachlich auf das Kapitel iiber die Betriebsparametertheorie, 
wobei die Abschnitte 13 bis 16 neu hinzukamen. Ebenfalls aus dem Nachlasse von 
CavER stammt der Abschnitt iiber Filter mit Elektronenréhren (VI./19, S. 275). 
Ein Einfiihrungskapitel, das den Stoff des ganzen Buches zusammenfassend 
behandelt, erleichtert dem Leser die Ubersicht. Die Anhange zu den einzelnen 
Kapiteln in der 1. Auflage wurden an den Schluss des Buches verwiesen, ebenfalls 
die fiir den Ingenieur sehr wertvollen Hilfsmittel fiir den Entwurf von Filtern, ge- 
trennt nach der Wellenparametertheorie und nach der Betriebsparametertheorie. 
Die Erlauterungen zum Entwurf von Filtern nach der Wellenparametertheorie 
wurden von den Herausgebern iiberarbeitet und erweitert, so dass es auch den 
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mit der Theorie weniger Vertrauten mit diesen Hilfsmitteln mOglich sein sollte, 
Filter zu berechnen, die bestimmten Anforderungen geniigen. Fiir Leser, welche 
das Buch von CavER noch nicht kennen, folgen einige Titel von Kapiteln: 1. Auf- 
gabenstellungen und Beispiele, 2. Berechnung der Eigenschaften gegebener Schal- 


tungen (unter anderem Topologie), 3. Vierpole, 4. Positive Funktionen und positive 


Matrizen, 5. Reaktanztheoreme, 6.Wellenparametertheorie der Tiefpass-Reaktanz- 
filter, 7. Wellenparametertheorie allgemeiner Filter, 8. Reaktanzvierpole mit vor- 
gegebenen Betriebseigenschaften, 9. Frequenzweichen, 10. Aquivalenz von Re- 
aktanzschaltungen. Anhange: I. Hilfsmittel aus der linearen Algebra, II. Elemente 
der Theorie der analytischen Funktionen, III. Vierpol-Formelsammlung, 1V. L6- 
sung einiger Tschebyscheffischer Extremalprobleme, V. Hilfsmittel fur den prak- 
tischen Entwurf von Filtern nach der Wellenparametertheorie, VI. desgleichen von 
Filtern und Frequenzweichen nach der Betriebsparametertheorie. 

Druck und Einband entsprechen dem innern Wert des Buches, sind also gut. 
Druckfehler sind bei einer solchen Fiille des Stoffes kaum zu vermeiden und 
darum entschuldbar. Mége diese Neuauflage eine weite Verbreitung finden als 
Unterstiitzung der Ingenieure in der Netzwerk-Analyse und Synthese. Zugleich 
ist das Buch ein Beweis dafiir, was angewandte Mathematik zu leisten vermag. 

H. Weber 


Existence Theorems for Ordinary Differential Equations. Von F. J. Mur- 
RAY und K. S. MirteR (New York University Press, New York 1954). 164 S., 
6 Fig.; $5.00. 

Das Buch umfasst die klassischen Theoreme und bringt zunachst den funda- 
mentalen Existenzsatz von Peano und den Eindeutigkeitssatz fiir Systeme von 
beliebig vielen gew6hnlichen Differentialgleichungen in expliziter oder impliziter 


Form; fiir dieselben Falle wird mit Hilfe der Methode der sukzessiven Approxi-. 


mationen die Stetigkeit in den Anfangswerten und Parametern nachgewiesen. 
Alsdann werden die Eigenschaften der Losungen, wie Differenzierbarkeit nach den 
Anfangsbedingungen und Parametern, untersucht und schliesslich kurz die 
linearen Differentialgleichungen im Reellen und Komplexen betrachtet. Die Dar- 
stellung ist tibersichtlich und klar, die Voraussetzungen werden sorgfaltig hervor- 
gehoben und die Beweise auch im allgemeinen Falle ausfiihrlich dargestellt. Die 
Verfasser weisen auf die Bedeutung dieser Satze bei der Lésung von Differential- 
gleichungen mit Hilfe von Rechenmaschinen (Analogiegeraten) hin. — Das Buch ist 
in erster Linie fiir Studenten in mittleren Semestern gedacht. E. Roth-Desmeules 


Die Messwandler, ihre Theorie und Praxis. Von J. GoLpsTEIN (Verlag 
Birkhauser, Basel 1952). 2. Auflage, 224 S., 210 Fig.; Ganzleinen Fr. 29.10, bro- 
schiert Fr. 24.95. 

Der Autor war friiher bei der AEG., Berlin, und spater als beratender Inge- 
nieur der Firma Moser-Glaser & Co. in Muttenz als anerkannter Spezialist auf dem 
Wandlergebiet tatig. Die hier vorliegende 2. Auflage beriicksichtigt die Weiter- 
entwicklung der Wandler und behandelt auch die amerikanischen Arbeiten. 

Zuerst wird die Theorie des Stromwandlers, der Stromfehler und der Fehl- 
winkel an Hand von Zeigerdiagrammen klargestellt, anschliessend kommen die 
Materialeigenschaften der modernen magnetischen Blechsorten, die Dimensionie- 
rung der Kerne, die Berechnung und Messung der Streuung zur Behandlung. Die 
nachsten Kapitel beschaftigen sich mit den Mitteln zur Verbesserung der Wandler- 
eigenschaften, namlich mit dem Zweistufenwandler, der Wilson-Kompensation, 
der Methode von A.C. ScHwacer, dem kompensierten und insbesondere dem 


a 
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vom Autor geférderten vormagnetisierten Wandler. Es folgen nun die fiir die 
Praxis wichtigen Probleme des konstruktiven Aufbaues, der dynamischen und 
thermischen Kurzschlussfestigkeit, der Uberstromziffer und der StoBspannungs- 
festigkeit. Der Gleichstrommesswandler wird nur kurz gestreift, da fiir diesen 
noch keine Fehlertheorie vorliegt. Den Abschluss bilden Betrachtungen iiber den 
Spannungswandler und die Prazisionsmesstechnik fiir Wandlerpriifungen. 

Der Text ist klar geschrieben und durch viele Zeigerdiagramme, Kurven und 
Photographien erlautert. Da das Buch sich an einen grésseren Kreis von Elektro- 
ingenieuren, Messtechnikern und Betriebsingenieuren richtet, ist von der kom- 
plexen Rechnung kein Gebrauch gemacht worden, obschon dadurch gewisse 
Uberlegungen gewonnen hatten. Die drucktechnische Ausstattung ist hervor- 
ragend. Leider war bei der Drucklegung die Publikation Nr. 192df des SEV. noch 
nicht erschienen, so dass von den SEV.-Symbolen und ferner vom Giorgi-System 
noch kein Gebrauch gemacht werden konnte. Das Buch kann warm empfohlen 
werden. Ed. Gerecke 


Siebenstellige Tafeln der elementaren transzendenten Funktionen. 
Von Dr. FriepRicH Léscu (Springer-Verlag, Berlin 1954). 331 S.; DM. 49.80. 

Zu jedem Argument x (Bogenmass) sind tabelliert die Funktionen arc x, 
sin %, cos ¥, tg ¥, arcsin x, arccos #, ln x, e”, e-® und die hyperbolischen Funktionen 
sinhw, cosh, tgh#, arsinh#, artgh#, und zwar in Tafel I fiir ¥ = 0 (0,0001) 0,1 
neunstellig und in Tafel II fiir ¥ = 0,1 (0,0005) 3,15, ¥ = 3 (0,01) 10 und ¥ = 10 
(0,1) 20 siebenstellig. Die ersten Differenzen sind durchweg mitgedruckt, wobei 
in konsequenter Weise als Schritt immer eine Einheit der letzten Stelle von * 
zugrunde gelegt ist. In allen Tafeln ergibt lineare Interpolation einen Fehler 
von héchstens 2 Einheiten in der letzten Stelle des Funktionswerts. Verschiedene 
weitmaschigere Hilfstabellen beschliessen das Buch, das wegen seiner tibersicht- 
lichen Anordnung und bequemen Handhabung verdient, als Standardwerk der 
elementaren Transzendenten zu gelten. E. Stiefel 


Dictionary of Conformal Representations. Von H. KoBER (Dover Publi- 
cations 1952) 208 S., 447 Fig.; $3.95. 

Dieses sehr niitzliche Handbuch erfiillt einen doppelten Zweck. Es enthalt 
einerseits in zeichnerischer Darstellung die konformen Abbildungen, welche durch 
gegebene analytische Funktionen von einfachem Typus hergestellt werden, und 
es erlaubt andererseits, die konforme Abbildung von gegebenen Bereichen auf- 
einander zu konstruieren, insofern dies durch geschlossene Ausdriicke méglich 
ist. Behandelt werden: Die ganze und die gebrochene lineare Funktion, quadra- 
tische und Potenzfunktionen, Exponentialfunktion und Logarithmus, elliptische 
und Modulfunktionen. Hervorgehoben sei die ausfiihrliche Behandlung der durch 
elementare Funktionen realisierbaren Schwarz-Christoffelschen Abbildungen 
sowie die Beispiele, welche die Abbildung von Tragfliigelprofilen betreffen. Das 
Buch ist fiir die Behandlung von Randwertaufgaben der mathematischen Physik 
und speziell fiir die Stromungslehre ein ausgezeichnetes und zeitsparendes Hilfs- 
mittel. E. Stiefel 


Numerical Solution of Differential Equations. Von W. E. MILNE (John 
Wiley & Sons, New York 1953) 275 S.; $6.50. 

Wer oft mit Differentialgleichungen zu tun hat, weiss, wie haufig man auf 
numerische Lésungsmethoden angewiesen ist. Das vorliegende Lehrbuch gibt 
einen Querschnitt durch das gesamte Gebiet der numerischen Behandlung von 
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Differentialgleichungen. Dabei ist es dem Verfasser gelungen, den Umfang zu © 


beschranken, ohne auf die Behandlung der zahlreichen rechentechnischen Kunst- 
eriffe verzichten zu mussen. ' 
Nach einem einleitenden Kapitel iiber die Grundbegriffe der numerischen 
Integration behandelt der Verfasser im 2. Kapitel die einfachsten Integrations- 
methoden (zum Beispiel die Trapezregel). Es folgt ein Kapitel 3 tiber theoretische 
Grundlagen, worauf in den Kapiteln 4, 5, 6 die verfeinerten Methoden zur nume- 
rischen Integration von Differentialgleichungen und Differentialgleichungssy- 
stemen zur Sprache kommen. Kapitel 7 ist den Randwertproblemen gewidmet, 
wobei auch die Verfahren von Ritz und GALERKIN gestreift werden. Der zweite 
Teil des Lehrbuches befasst sich mit partiellen Differentialgleichungen, mit einer 
starken Betonung der elliptischen Gleichungen. Demgemass ist auch ein betracht- 
licher Teil (Kapitel 9) den weiteren Ausfiihrungen der Losung linearer Gleichungs- 
systeme gewidmet. Auch die Bestimmung von Eigenwerten und Eigenvektoren 
einer Matrix kommt in Kapitel 9 zur Sprache, obwohl sich das Lehrbuch sonst 
nur wenig (in Kapitel 11) mit Eigenwertproblemen befasst. Das Werk schliesst 
mit einem iiberaus reichen Literaturverzeichnis. HA, Rutishauser 


Handbook of Elliptic Integrals for Engineers and Physicists. Von 
PauL F, Byrp und Morris D. FRIEDMAN (Springer-Verlag, Berlin 1954). 355 S., 
22 Fig.; DM. 36.—/39.60. . . 

Den Hauptinhalt des Buches bilden etwa 3000 Formeln zur Berechnung ellip- 
tischer Integrale, wobei sehr breit die verschiedenen Spezialfalle des unter der 
Wurzel stehenden Polynoms vierten Grades beriicksichtigt sind. Das Buch leistet 
daher niitzliche Dienste zur Abkiirzung der an sich elementaren Reduktion eines 
elliptischen Integrals auf Normalform. Die Verfasser benutzen ausschliesslich die 


Jacobischen elliptischen Funktionen fiir ihre Formeln; zur Berechnung der — 


Integrale dritter Gattung wird ausserdem die Heumansche Funktion herange- 
zogen, die zusammen mit den Normalintegralen erster und zweiter Gattung 
- sowie der Z-Funktion am Schluss des Buches kurz tabelliert ist. Die WeierstraB- 
schen Funktionen und die #-Reihen werden nur im Anhang kurz beriihrt, sind 
aber nicht in die Formelsammlung eingebaut. Das Buch ist so angelegt, dass es 
ohne tieferes Eindringen in die Theorie der elliptischen Funktionen benutzt 
werden kann. E, Stiefel 


Projektive Geometrie. Von W. BLAscHKE (Birkhauser Verlag, Basel 1954), 
197 S. mit 71 Abb., geb. sFr. 19.60, brosch. sFr. 16.65. 

Die Kenntnis der Figuren und Satze der projektiven Geometrie ist unerliss- 
liche Voraussetzung fiir jede Beschaftigung mit irgendeinem Zweige der Geome- 
trie. Die projektive Geometrie gehért daher zum Unterrichtsstoff der unteren Se- 
mester des Mathematikstudiums. Aus einer derartigen Vorlesung ist BLASCHKES 
«Projektive Geometrie» hervorgegangen. 

Zur Behandlung der projektiven Geometrie stehen zwei Wege offen, der ana- 
lytische und der synthetische. Dem propddeutischen Charakter des Buches ent- 
sprechend wahlt der Verfasser die analytische Methode. Als Koordinatenk6érper 
ist der K6rper der komplexen Zahlen zugrunde gelegt. 

Gegenstand der projektiven Geometrie ist die Untersuchung jener Eigen- 
schaften von Figuren, die gegentiber Kollineationen invariant sind. Dass bei 
einer Einschrankung der Transformationsgruppe der Invariantenschatz grésser 
wird, fihrt der Verfasser an einigen Beispielen ausfiihrlich vor. Insbesondere 
werden die Affingeometrie, die Euklidische Geometrie und die Nichteuklidische 
Geometrie unter diesen Gesichtspunkten behandelt. 
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Verschiedentlich werden auch Briicken zur darstellenden Geometrie geschla- 
gen, indem etwa ebene Figuren als Projektionen gewisser raumlicher Figuien 
interpretiert werden. Der Leser gewinnt dadurch auch einen kleinen Einblick in 
die synthetische Methode. Weitere Querverbindungen reichen zur algebraischen 
Geometrie, zur graphischen Statik und zur Algebra. Dazu gehért insbesondere die 
Aufnahme des von-Staudtschen Hauptsatzes iiber die Abbildungen, welche die 
harmonische Lage von 4 Punkten bewahren. 

Das Buch gliedert sich in folgende Abschnitte: I, Einleitung, II. Homogene 
Koordinaten, Kollineationen und Korrelationen, III. Das Doppelverhaltnis 
IV. Kegelschnitte, V. Liniengeometrie, VI. Quadriken, VII. Nichteuklidische 
Geometrie, VIII. Mébiussche Vierflachpaare, IX. Gewebe und Gruppen. 

Im letzten Teil tiber Gewebe und Gruppen werden einige axiomatische Fragen 
der projektiven Geometrie gestreift unter Verwendung neuerer Begriffsbildungen. 

Die «Projektive Geometrie» BLASCHKEs ist 4usserst ansprechend geschrieben 
und kann zur Einfiihrung in dieses Gebiet bestens empfohlen werden. Neben dem 
Neuling wird aber auch der Kenner mit Vergniigen darnach greifen, enthalt doch 
das Buch eine Reihe von weniger bekannten geometrischen Zusammenhangen. 

M. Jeger 


Nutshell Tables of Mathematical Functions for Interpolation with 
Calculating Machines. Von R. A. HiRVoNEN. 24 S. Sonderdruck aus Bull. géo- 
désique 30 (1953). 

Auf 4 Druckseiten zusammengedrangte 10stellige Tafeln der sogenannten ele- 
mentaren Funktionen (log x, 10”, cos %, sinx, arctg x) sowie 3 Seiten Tabellen fiir 
die Berechnung geodatischer Linien. 

Der Verfasser teilt das Argument in relativ grosse Intervalle ein und gibt dann 
fiir jedes dieser Intervalle ein Polynom vierten Grades an, welches dort um 
héchstens 0,5-107-19 von der betreffenden Funktion abweicht. Das Aufschlagen 
eines Funktionswertes ist also mit der Auswertung eines Polynoms vierten 
Grades verbunden, was aber bei Vorhandensein einer kleinen Rechenmaschine 
durchaus ertraglich ist. H. Rutishauser 


Mathematics and Plausible Reasoning. Von GreorceE Potya. I. Induction 
and Analogy in Mathematics. 11. Patterns of Plausible Inference (Princeton Uni- 
versity Press, Princeton, New Jersey, 1954). Bd. 1: 292°S., 64 Fig.; $5.50. 
Bd. 2: 200 S., 13 Fig.; $4.50. 

Es gehdért seit den Zeiten EuKiips zum Charakter der Mathematik, dass die 
Darstellung eines Sachverhaltes in eine méglichst unpers6nliche und in gewissem 
Sinne endgiiltige Form gebracht wird. Mit der gréssten Sorgfalt werden die 
Spuren, die zu einem Satz oder seinem Beweis fiihrten, ebenso alle Umwege ver- 
wischt, und vieles erscheint dann als ein «Deus ex machina». Es gibt nur wenige 
Mathematiker, man denke etwa an EvLeEr, die in aller Offenheit auch ihre un- 
fruchtbaren Versuche erwihnen. Dass gerade das Finden des Weges zur Losung 
eines Problems das Wichtigste ist, ist offensichtlich; aber diese Seite der Mathe- 
-matik wird in den allermeisten Biichern iiber Mathematik vernachlissigt. 

Nun hat G. Potya seine Erfahrungen auf dem Gebiete des plausiblen Den- 
kens, iiber das er schon an den verschiedensten Orten Vorlesungen hielt und eine 
Reihe von Abhandlungen ver6éffentlicht hat, in einem ausgezeichneten, zweiban- 
digen Werk zusammengefasst. Das erste Buch zeigt an Hand sehr instruktiver 
Beispiele aus der elementaren und héhern Mathematik und der mathematischen 
Physik die Bedeutung auf, welche die Analogie und die Induktion im allgemeinen 
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Sinne (das induktive Denken spielt ja in den Naturwissenschaften und im tag- 
lichen Leben eine hervorragende Rolle) beim Auffinden von Vermutungen und 
schliesslich beim deduktiven Beweisen besitzen. Damit wird gleichzeitig der Grund 
gelegt fiir die Ausfiihrungen des zweiten Bandes. Dieser bringt zunachst eine 
Systematisierung der plausiblen Schlisse, ihre gegenseitigen Beziehungen, und 
nach einem Uberblick iiber die Wahrscheinlichkeitsrechnung folgt der bemer- 
kenswerte Versuch einer Verbindung derselben mit dem plausiblen Schliessen, 
indem an Stelle der Wahrscheinlichkeit die Glaubwiirdigkeit gesetzt wird. Selbst- 
verstandlich kann es sich dabei nur um qualitative Aussagen handeln, doch 
erscheint der Weg vielversprechend. — Allen Kapiteln sind sorgfaltigst ausge- 
wahlte Ubungen und Erganzungen zum Text beigefiigt. 

Das Werk von G. Potay, das eine Hilfe (a guide to the art of plausible 
reasoning) fiir Studenten und Lehrer sein will, ist gleichzeitig in einem gewissen 
Sinne ein philosophischer Essay und stellt einen wichtigen Beitrag zu dem wenig 
beachteten Gebiet der Psychologie des Erfindens dar. Es besteht kein Zweifel, 
dass das feine Werk, welches die jahrzehntelangen Erfahrungen eines gewiegten. 
Mathematikers, Problemstellers und -lésers enthalt, vor allem den Studenten, 
aber auch dem Praktiker in seinen Mussestunden tiefe und wertvolle Einsichten. 
vermittelt. So ist zu hoffen, dass das Buch die weiteste Verbreitung finden und 
die Bemiihungen des Verfassers reiche Friichte tragen werden. 

E. Roth-Desmeules 


The Statistical Approach to X-ray Structure Analysis. Von V. VAND- 


und R. Pepinsky (The Pennsylvania State University 1953). 109 S., 18 Fig.; © 


$1.50. 
Die Monographie ist im wesentlichen eine kritische Betrachtung des Werkes. 


von H. HAUPTMANN und J. Kare: Solution of the Phase Problem I. The Centro- - 


symmetric Crystal, ACA Report 3 (The Letter Press, Wilmington, Del., 1953), in 
welchem eine statistische Behandlung des Phasenproblems in der Kristallstruk- 
turanalyse unternommen worden war. Die Autoren zeigen die Korrelation der 


abgeleiteten statistischen Formeln mit Modifikationen von Patterson- und Pat-— 


terson-Harker-Analysen und den Vorzeichenkriterien von SAYRE-COCHRAN- 
PATTERSON und weisen am Beispiel einfacher Punktstrukturen und durch all- 
gemeine Uberlegung nach, dass die genannten statistischen Formeln keine all- 
gemeine Giiltigkeit besitzen k6nnen. W. Hoppe 


Vorlesungen tiber Baustatik, Vol. IJ. Par Fritz Srutssi (Birkhauser, 
Basel 1954). 313 p., 217 fig.; reliure toile fr. 37.45, broché fr. 33.30. 

Dans le premier volume de son «Cours de statique des constructions», dont la 
premiére édition est parue en 1946 et la seconde en 1953, le professeur FRITz. 
STUssI expose les théories et les méthodes permettant de calculer les tensions et 
les déformations des systémes statiquement déterminés. 

Le second volume, qui vient de paraitre et que nous présentons aujourd’hui. 


aux lecteurs du ZAMP, traite avant tout les systémes statiquement indéterminés. — 
Il s’agit donc d’un ouvrage qui vient heureusement compléter la matiére du 


tome I. 

Le premier chapitve de ce nouveau livre est consacré aux propriétés fondamen- 
tales des déformations élastiques des éléments de construction situés dans un 
plan. Dans le second chapitre, auteur expose les principales méthodes permettant 
de calculer les systémes statiquement indéterminés. Avec juste raison, il insiste 


sur la méthode classique consistant a remplacer le systéme étudié par un systéme 
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fondamental convenablement choisi, sur lequel sont appliquées non seulement 
les forces données, mais encore les grandeurs surabondantes. Tl décrit aussi cn 
détail les méthodes de calcul permettant de résoudre élégamment, par rapport a 
ces grandeurs, les équations régissant les déformations. 

Les autres chapitres (III & VII) sont consacrés a l’application a divers systemes 
statiquement indéterminés, des méthodes décrites dans la premiére partie. I] 
S’agit principalement de la poutre continue sur appuis fixes ou élastiquement 
déformables, des avcs a deux et a trois articulations, de Varc encastvé, de différents 
types de cadres, y compris la poutre Vierendeel, des poutres venforcées par une 
chainette et des systémes triangulés statiquement indéterminés. Quelques considéra- 
tions sur les tensions secondaires dues a la rigidité des attaches des barres dans les 
systémes triangulés terminent l’ouvrage. 

Comme le premier volume, ce livre est trés bien présenté. I] est écrit dans un 
style clair, précis et agréable a lire, I] rendra donc service non seulement aux 
étudiants des écoles techniques de degré universitaire, mais aussi aux ingénieurs 
de la pratique, qui seront heureux d’y trouver d’abord une remarquable vue 
d’ensemble sur les questions traitées, puis un outil précieux pour l’élaboration de 
leurs projets. H, Favre 


Theoretische Bodenmechanik. Von K.TERZzAGHI, iibersetzt von R. JELINEK. 
(Springer-Verlag, Berlin 1954). 505 S., 153 Abb.; Ganzleinen DM. 48.-. 

Das Buch iiber «Theoretical Soil Mechanics» erschien 1942 in den USA. Es 
fand bald Eingang in weiteste Fachkreise und stellt heute eines der massgebenden 
Werke in der Erdbaumechanik dar, besonders auch deshalb, weil es fiir den Bau- 
ingenieur geschrieben worden ist. Die Ubersetzung von Dr. JELINEK schliesst nun 
auch im deutschen Sprachgebiet eine oft empfundene Liicke. 

Die «Theoretische Erdbaumechanik» fiihrt den Leser nicht, wie auf Grund 
des Titels zu vermuten ware, in das rein abstrakte Wissensgebiet der Erdbau- 
mechanik, sondern sucht vielmehr eine Briicke zwischen Theorie und Praxis zu 
schlagen. Der Verfasser weist darauf hin, dass in Analogie zu den andern Fach- 
gebieten der Ingenieurwissenschaften die bodenphysikalischen Eigenschaften 
weitgehend idealisiert worden sind, so dass deren Verhalten in Form einfacher 
Gesetze wiedergegeben werden kann. Um die Ubersichtlichkeit zu wahren, 
wurde auf die Wiedergabe von Herleitungen verzichtet und manche numerisch 
nicht lésbare Probleme durch Naherungen ersetzt. 

Das Buch behandelt folgende Fachgebiete: 

1. Allgemeine Grundlagen, umfassend die Grundelemente, die plastischen 
Grenzzustande im isotropen Halbraum und theoretische Grundlagen fiir die 

-Lésung praktischer Aufgaben. Die Mohrsche Bruchhypothese findet weitgehende 
Anwendung. 

2. Brucherscheinungen in idealen Béden, Erddruck und Erdwiderstand, Trag- — 
fahigkeit, Standsicherheit von Béschungen und Spundwande. Von besonderer 
Bedeutung fiir das Bauwesen ist das Kapitel iiber die Tragfahigkeit der Béden. 
Seine Anschauungen waren richtunggebend fiir die Entwicklung der Fundations- 
technik. Es ist dem Verfasser hoch anzurechnen, dass er aus den vielenTheorien 
besonders iiber das Gebiet der Standsicherheit von Béschungen in klarer Form 
die massgebenden Faktoren herausgearbeitet hat, um so dem weniger Fach- 

-kundigen einen guten Uberblick zu verschaffen. 

3, Mechanische Wirkung des Wassers im Boden. Behandelt werden die Sicker- 
str6mung und deren Einfluss auf die Beanspruchung des Bodens, Setzungsvor- 
gange (Konsolidationstheorie) und Kapillarerscheinungen. Die von K. TERZAGHI 
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und O. K. FROHLICH ausgearbeitete Konsolidationstheorie eonde besondere , 
Beachtung, da es heute méglich ist, damit die Spanniungsverhiltnisse im Bau-_ 
grund wahrend eines Setzungsvorganges abzuschatzen,. 

4, Elastizitatsaufgaben der Erdbaumechanik unter Beriicksichtigung der - 
Spannungsverteilung im elastisch-isotropen Halbraum, der Theorie elastischer 
Schichten auf starrer Unterlage und dynamischer Lasten. 

Obwohl seit dem Erscheinen der 1. Auflage (1942) die Erkenntnisse enorme 
Fortschritte gemacht haben, verlor das Buch keineswegs an Wert, wurde es ~ 
doch fiir die folgenden Auflagen laufend erganzt und dem neuesten Stand ange- 
passt. Der besondere Vorzug dieses Buches liegt in der leichten Uberblickbarkeit 
der Materie und der kritischen Betrachtung der Theorien. 

Die Ubersetzung von Dr. JELINEK ist prazis und halt sich genau an den- 
amerikanischen Originaltext. Es ist nur zu hoffen, dass auch das 2. Buch von 
Prof. TERZAGHI: Soil Mechanics in Engineering Practice, im deutschen Sprach- 
gebiet Eingang findet. i ¥eZellers 


Electron Optics. Von O. KLEMPERER, 2. Aufl. (Cambridge University Press 


“' 4953)..471 S., 167 Fig., 50 s. 


». Diekerste ‘Auflage ist 1939 erschienen. Dass der Umfang von etwa 100 Seiten. 
dieser ersten Auflage auf iiber 450 Seiten der zweiten Auflage angewachsen ist, i 


- wird nicht iiberraschen. Einmal, weil die technische Bedeutung der Elektronen-— 


optik erheblich grésser geworden ist; dann aber auch, weil in den seit der ersten © 
Auflage verstrichenen Jahren sehr viele neue und wesentliche Erkenntnisse rein 


-. wissenschaftlicher Art hinzugekommen sind. 


Es ist deshalb eher erstaunlich, dass es dem Verfasser setimzen ist, ‘aed dem 


-relativ geringen Raum von 450 Seiten eine Darstellung zu geben, welche ein recht 
‘vielseitiges Bild vom heutigen Stand der experimentellen und theoretischen 


Elektronenoptik zu vermitteln vermag. Allerdings ist das Schlusskapitel iiber 
Anwendungen der Elektronenoptik summarisch geraten und eher als ein erster 
Hinweis zu betrachten. Demgegeniiber sind alle anderen Kapitel eine stellenweise 
sehr ins einzelne gehende Beschreibung, die der zweifachen Zielsetzung des — 
Buches ~ es soll einerseits eine Einfiihrung fiir Studierende sein, andererseits dem — 


- Fachmann die numerischen Unterlagen fiir seine Entwicklungsarbeit vermitteln — 


voll gerecht wird. Der Schwerpunkt der Darstellung liegt mehr auf der experi- 
mentellen Seite. Auf die Wiedergabe der Zwischenrechnungen wird haufig ver- 
zichtet. Es mag vielleicht als ein Nachteil bewertet werden, wenn der Leser 


4 gewissen Formeln gegeniibergestellt wird, ohne in der Lage za sein, den Rech- 
nungsgang zu wiederholen; es geht damit der Darstellung eine gewisse Zwangs- 


laufigkeit verloren. Jedoch nur auf diese Weise war es moglich, den Umfang des 
Buches derart zu halten, dass der Leser den Uberblick iiber das Wesentliche nicht | 


nt verliert. Jedenfalls kennen wir heute kein anderes neueres Werk, welches_ als 


erste Einfiihrung in die experimentelle Elektronenoptik besser géeignet ware. 

Das Buch umfasst folgende Kapitel: Historical Introduction: The Fundamen- 
tal Principles of Electron Optics. The Cardinal Points of an Electron Lens. Field 
Plotting and Ray Tracing. Some Electrostatic Electron Lenses. Magnetic Electron 
Lenses. Lens Errors: Geometrical Aberrations. Electronic Aberrations. Electron 
Optics and Space Charge. Emission Systems. Electron Lenses and. Emission 
Systems with Line Focus. Deflecting Fields. Application of Electron Optics in 


? Industry and Reséarch. Ein Verzeichnis von etwa 500 Verdffentlichungen sowie 


eine Anzahl von Tabellen sind beigefiigt. A, o: Rusterholag 
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